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Abstract 

A simple, rapid, sensitive and accurate spectrophotometric method is developed for the 
determination of carbosulfan and propoxur in field water samples and grains The method is 
based on the formation of coloured species upon coupling their hydrolysis products with 
diazotized p-aminoacetanilide, exhibiting maximum absorption at 390 and 389 nm respectively. 

The method obeys Beer’s law over the concentration range 0 5-12 and 0.5-16 pg/ml. 

(Keywords, spectrophotometry/carbosulfan, propoxur determination / p-amino-acetanilide / 
chloroform) 


Introduction 

Carbosulfan, carbamic acid [(dibutyIamino)thio]methyl-2,3-dihydro-2,2-dimethyl- 
7-benzofuranyl ester is a recently developed broad spectrum systemic and contact 
insecticide belongs to the class of carbamate insecticides. It is used for the control of 
soil dwelling insects and foliar pests on maize, potatoes and sugarbeet and also used 
as a nematicide. It acts as a cholinesterase inhibitor. Carbosulfan is unstable in acid 
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medium and it will be converted into carbofuran by N-S bond cleavage. Clay et al\, 
have reported the activity of this compound 

Propoxur, 2-(l -methyl ethoxy) phenyl methyl carbamate also belongs to 
carbamate class of pesticides. It is a white crystalline solid having m.p. 91”C, soluble 
in most organic solvents . Propoxur is commercially available as 1% Baygon spray, 
2% bait, 4% Baygon and 70% wettable powder. Baygon has rapid knockdown and 
long residual action. It shows cholinesterase inhibiting properties and is particularly 
effective as an insecticide. It is particularly effective against insects affecting man and 
animals, such as cockroaches, flies and mosquitoes. Foilage sprays have given 
effective control of aphids, lygus bugs, grass hoppers and other insects on various 
crops. Bioactivity and persistence of propoxur in mineral soil, river water and in the 
aquatic environment have been reported^''*. Wide spread application of these 
carbamate pesticides to control pests in agricultural as well as domestic applications. 
All these pesticides from the agricultural field, and domestic applications are washed 
over the nearby water sources either by leaching or by the surface runoff. Such residue 
concentration are more than the permissible levels and are generally harmful to pond 
fauna and the people who drink these waters. 

Literature survey revealed that very few analytical methods such as GC , 
HPLC and TLC and some other spectrophotometric methods’^'^^ have already 
been reported. The present communication reports a new spectrophotometric method 
for the determination of the above pesticides based on the formation of a coloured 
species by coupling the corresponding phenols obtained by alkaline hydrolysis with 
diazotized p-aminocetanilide. 


Materials and Method 

All chemicals used were of analytical grade and all of the solutions were freshly 
prepared with distilled water. A Hitachi, model U-3400, UV-VIS NIR 
Spectrophotometer with 1 cm glass cells was used for absorbance measurement. 
Sodiumhydroxide (2%) and sodmmnitrite (0.5%) solutions were freshly prepared, p- 
aminoacetanilide solution (0.1%) is prepared in 100ml IN hydrochloric acid. 

Preparation of standard solutions of Carbosulfan and Propoxur (50fj.g/ml) : 

Stock solution of the above pesticides were prepared by dissolving 50 mg of 
analytical grade insecticide in 100ml of methanol. From this working standard 
solution is prepared by diluting 10 ml of this solution to 100ml. 
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Preparation of sample solution . 

Formulations; Carbosulfan 25% and Propoxur 70% wettable powder. An amount 
equivalent to 50 mg of the insecticide was dissolved in 50 ml of methanol. 5 ml 
aliquots of the solution was subsequently diluted with 100ml of methanol. 

General Procedure : 

Aliquots of standard insecticide solution (50|ig/ml) 0, 0.25, 0 5, 1.0, , 6m! 

were placed in series of 25 ml standard flasks followed by adding 2.5ml of 2% 
sodiumhydroxide and 2.4ml of diazotization mixture. The solutions were made to the 
mark with distilled water and absorbance was measured for the yellow coloured 
chromophore against a reagent blank. The absorbance is plotted against the final 
concentration to obtain a calibration graph. 

Water samples ■ pH of each water sample was adjusted to 3-4 with 20% sulphuric 
acid. One litre samples of distilled water and tap water was fortified with different 
concentrations of insecticide dissolved in methanol The fortified water samples were 
extracted with 100ml chloroform using a seperating funnel. Chloroform extracts were 
collected into a funnel and re-extracted the aqueous phase twice with 50ml 
chloroform. The combined extracts were washed with 0.1 M potassium carbonate 
solution and dried over anhydrous sodium sulphate in a filter funnel and collected the 
extracts in a 250ml standard flask. Finally these were made upto the mark with 
chloroform. Known aliquots of the chloroform were evaporated to dryness on a steam 
bath. The residue was dissolved in methanol, then developed the colour. 

Grains : lOOg of grains (rice and wheat) was taken in a conical flask and shaken 
for 5 min, with 200ml chloroform. Chloroform was filtered into a 250ml standard 
flask through a Whatman No.l filter paper and residue was washed twice with 10ml 
chloroform. Chloroform extracts were combined and made up to the mark. Known 
aliquots of the chloroform extracts were used for colour developement after 
evaporating chloroform. 

Results and Discussiotii 

The absorbance maximum was at 3S!0nm for Carbosulfan and 389nm for 
Propoxur. Beer's law is obeyed over the range 0.5-12 pg/ml for Carbosulfan and 
0.5-16 |J.g/ml for Propoxur. The colour developes instantenously and remains stable 
for more than 24 h. The coupling reaction can be done at room temperature only. 

The suitability of the proposed method was studied by analysis of ten replicate 
samples containing 5ppm of Carbosulfan and Propoxur The relative error and relative 
standard deviations are given in Table 1. 
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Carbosulfan with p-aminoacetanilide 



CARBOSULFAN & PROPOXUR IN ENVIRONMENTAL SAMPLES 

Reaction Mechanism 




NHC0CH3 NHCOC113 


p-aminoacet anilide 

Propoxur with p-aminoacetanilide 

Table 1- Characteristics of the method 


Parameter 

Carbosulfan 

Propoxur 

Present work 

Reference 

method*^ 

Present work 

Reference 

method*® 

^maxs 

390 

475 

389 

500 

Beer's law limit, lig/ml 

0.5-12 

1-10 

0 5-16 

0 25-5 0 

Sandell's sensitivity, |ig cm'^ 

0 0083 

- 

0 0125 

- 

Molar absorptivity, litre mole‘^ cm"* 

4 576x10^ 

- 

1 674a 10^ 

- 

Stability of the coloured species, hr 

24 

8 

24 

12 

Relative Standard Deviation, % 

0.91 

0.8 

04 

06 

Relative Error, % 

0.5 

06 

02 

0.3 
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Formulations containing Carbosulfan and Propoxur were analysed (10 replicates). 
For a 25% Carbosulfan emulsion the mean ± standard deviation was 24.48 ± 0.12%. 
For a 1% Propoxur spray and 4% dust the corresponding values were 0.96 ± 0 08% 
and 3 84 ± 0.02% respectively. 

Recovery experiments were performed with known amounts of the compounds 
added to different samples of grains and water. For water samples a methanol solution 
of the insecticides was added Grain samples were spiked by adding a methanol 
solution of the insecticide to the dry grains and evaporating the solvent. 

The results presented in Table 2 show that recovery was in the range 95 to 97%. 
The results in Table 3 suggest that the method is applicable to the analysis of field 
water samples also. 


Table 2- Recovery of Carbosulfan and Propoxur from grains and spiked water samples. 




Carbosulfan 



Propoxur 




Recovery, % 



Recovery, % 


Sample 

Added, Present work 

ppm 

Reference 

method'^ 

Added, 

ppm 

Present work 

Reference 

method^® 

Rice 

1 0 

96 4±1 9% 

96 01 1.6 

1 0 

96 5 1 1 6 

96 011 2 


30 

96 2±10% 

96 0 1 1 2 

30 

96.25 1 1 0 

96 0 1 1 0 


50 

95 4±09% 

94 8 1 0 8 

50 

96 1710 8 

95 8 1 0 8 


70 

94 8 ± 0 7% 

93 510 6 

70 

95.38106 

94 7 1 0 6 

Wheat 

10 

96 40 ±1 2% 

95 0 1 1 4 

1.0 

97 01 1.4 

96 0 ±1 4 


30 

96 20 ±1 0% 

94.111 0 

3.0 

96 50 1 0 9 

9561 1 0 


50 

95 17 ±0.8% 

93 2 1 0 7 

5.0 

96 33 1 0 5 

94.7 1 0 9 


70 

95.0010 6% 

92 810.5 

7.0 

95 63 1 0.4 

930107 

Water 

10 

96 201 1.0% 

96.0 1 1 6 

1 0 

97.5 1 1 2 

97.2 ±1 0 


30 

96 75 1 1 0% 

96011 2 

30 

97.75 1 0 9 

97 0 1 1 0 


50 

960109% 

95 210.8 

50 

96 5 1 0 7 

96.3 1 0.9 


70 

95.60 1 0 7% 

94 610.7 

70 

95 5 1 0.6 

95 1 1 0 7 
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Table 3- Determination of Carbosulfan and Propoxur in field water samples 



Carbosulfan found, ppm 

Propoxur found, ppm 

Sample Volume, ml 

Present work 

Reference 

method’^ 

Present work 

Reference 

method*^ 

250 

0 20 

0 20 

0 18 

0 14 

250 

0 30 

0 29 

0 28 

0 26 


The results obtained by the proposed method are found to be superior to the 
reference methods’®*’®. The method can be compared favourably with the reported 
methods’®"’®. As of its favourable characteristics as regards its stability, sensitivity 
and selectivity, the method can find applications to the field water samples also. 
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Abstract 

4,4'- Dihydroxy-3, 3'-dipropylbiphenyl-o-phenylenediainine (DDPBPD) and its chelates with 
Mn (111), Fe (III), Zr (IV), Th (IV), VO(II) and U02(I1) have been synthesised. These chelates 
have been characterized on the basis of elemental analysis, magnetic susceptibility, IR, 
electronic spectra and TGA analysis The I.R spectral studies indicate that the ligand acts as a 
dibasic tetradentate and coordinates through phenolic OH and azomethine nitrogen The 
dependence of the electrical conductivity on the temperature has been studied over the 
temperature range 300-520 K and the chelates are found to show semiconducting behaviour 
Thermogravimetric studies indicate the presence of water molecules in the Fe (III), Zr (IV) and 
Th (IV) chelates 

(Keywords IR/electronic spectra/polyschiffbase/polychelates/TGA/electncai conductivity) 

lotroduction 

Polymeric co-ordinating reagents are a novel type of substance possessing a 
combine action of the physical properties of high polymer and the chemical properties 
of the attached reagents’. The facile regenerability, higher stability and operational 
flexibility of the co-ordination polymers have created considerable research interest in 
the recent years^. The present paper describes the synthesis and characterization of 
chelates derived from 4,4'-dihydroxy-3-3'-dipropylbiphenyl-o-phenylenediamine. 
(DDPBPD). 


Materials and Method 

All the chemicals used were of Analar grade. DMF, ethanol and methanol were 
used after distillation. 4,4'-Dihydroxy-3,3’-dipropylbiphenyl (DDPB) was prepared by 
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Friedel-Crafts reaction. Manganese(III) acetate dehydrate was prepared^ by the 
oxidation of Mn (OAc )2 4 H 2 O. 

Preparation of ligand ■ 

The ligand DDPBPD was prepared by refluxing a mixture of DDPB and o- 
phenylenediamine in IT molar ratio in dry ethanol for 3 h. The resulting pale yellow 
solid was filtered and washed several times with hot water followed by ethanol and air 
dried. The product was insoluble in water and common organic solvents but soluble in 
hot DMF The product was crystallised from hot DMF The purity of the ligand was 
checked by tic. Yield (5.0 g, 68 %), m p 322 °C. 

Synthesis of polychelates 

A solution of DDPBPD (0.01 mol) was prepared in DMF (25 ml), a suspension of 
the metal salts [Mn (OAc )3 2 H 2 O, ferric chloride, vanadyl sulphate, zirconyl chloride 
and nitrates of Th and UO 2 ] m DMF-EtOH (1 : 4 v/v) (0.01 mol, in 25 ml) were added 
dropwise with constant stirring. A saturated solution of ammonium acetate was also 
then added with rapid stirring. On refluxing the mixture on a water bath for 6-7 hours,, 
solid were separated out, washed with hot water, DMF, EtOH and acetone and dried , 
over fused CaCl 2 . 

Carbon, hydrogen and nitrogen analyses of the ligand and polychelates were 
carried out on a Carlo Erba 1108 at microanalytical section Central Drug Research 
Institute, Lucknow. The metal content of the polychelates were determined 
volumetrically and gravimetrically^ (using a mixture of concentrated hydrochloric, 
nitric, per chloric and sulphuric acids ) IR spectra in KBr pellets were recorded on a 
Perkin-Elmer 842 IR spectrophotometer. Electronic spectra were recorded on a PAS-2 
spectrophotometer at BARC Mumbai. The magnetic susceptibility measurements 
were carried out at room temperature by the Gouy method using Hg[Co(NCS) 4 ] as a 
calibrant. The diamagnetic correction were taken from a standard source. 
Thermogravimetric analysis was performed on a Perkin-Elmer TGS-2 thermobalance 
with heating rate 10° min"' in air atmosphere at RSIC, Nagpur. For the electrical 
conductivity measurements a pellet of the polymer under study was prepared by 
pressing a finely powdered and dried sample under hydrostatic pressure at 20,000 Ibf 
in ^ for two min. A uniform thin layer of silver paste was applied to both sides of the 
pellet, providing electrical contacts. The average diameter and thickness of each pellet 
were measured. The pellet was firmly pressed between two circular metal discs 
functioning as electrodes. The other ends of the electrodes were passed through the 
pellet holder for connections. The entire assembly was placed in furnace. It was 
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heated at a rate of 1 ‘^'C min'. The resistance of the pellet was measured by Zentech 
702 A high resistance meter. 


Results and Discussion 

All the polychelates are coloured amorphous solids, insoluble in water and 
common organic solvents. The analysis of the compounds agrees with expected 1:1 
(metal: ligand) composition. The solids do not melt sharply but with decomposition. 
All the chelates are stable in air for extended periods of time. 

In the IR spectra of ligand medium broad band at 3237 cm"' due to 
intramolecularly hydrogen bonded OH stretch, disappeared in the spectra of all 
polychelates Indicating that the deprotonation of OH and its involvement in co- 
ordination^ The shift of C-0 from 1240 to 1280 cm'* supports bonding of ligand 
through the phenolic oxygen. The ligand band at 1638 cm'' due to central C=N 
stretch, shifted to lower frequency 1590 cm"' on chelation confirming the coordination 
through azomethine nitrogen atom to the metal ions. The bands at 970 and 900 cm in 
the spectra of VO(II) and UOaCII) polychelates were attributed to V=0 and 0=U=0 
stretches respectively. New bands observed in the region 450-650 cm are assigned to 
M-N and M-0 modes of vibration respectively^. 

The Zr(IV), Th(IV) and UOaCII) polychelates are found to be dimagnetic as 
expected. The electronic spectra of U02(II) polychelate exhibits charge transfer bands 
in the regions 32500-33200 and 35300-36100 cm"' assigned to n -» n and 7t -> n 
transitions respectively’’*. Low intensity band appearing around at 20000 cm"' shows 
the presence of symmetrical 0=U=0 stretching frequency^. The VO(II) polychelate 
has magnetic moment 1.79 B.M. and exhibits bands at 11764, 17182 and 25974 cm"' 
due to dxy -> dxz, dyz, dxy -» dx’-y’ and dxy -> dz’ transitions, respectively, towards 
square pyramidal stereochemistry*. The Fe(III) polychelate shows a magnetic 
moment value 5.97 B.M. correspondence to five unpaired electrons. The electronic 
spectra of this polychelate exhibits bands at 23255, 18939 and 17791 cm"', which are 
assigned to the transitions ®Aig -> ''Tig, ^Aig -> *T 2 g and ®Aig -> ''Aig , ‘*Eig , 
respectively'". From the electronic spectra and magnetic studies it is evident that 
Fe(III) polychelate possesses high spin octahedral configuration. The electronic 
spectrum of Zr (TV) polychelate exhibits single intense charge transfer band in the 
region 24860-28360 cm"' which is consistent with («-l) d° ns° electronic 
configuration. The Mn(III) polychelate shows magnetic moment value 4.95 B.M. 
expected for a high spin complex with four unpaired electrons. The electronic spectra 
of Mn(in) polychelate exhibits intense and sharp band around at 22700 cm'' and 
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another set of band in the region 12000-18520 cm'' are observed. The band around at 
22500 cm'^ is assigned to a charge transfer because of its high intensity. This 
assignment is in agreement with those of previous workers"’ . The splitting of d-d 
transition band in to a number of bands in the region 12000-18520 ern may be 
explained on the basis of a strong Jahn-Teller effect expected for a (f ion and 
lowering of symmetry from octahedral. 

The thermograms of these polychelates show two steps decomposition. The initial 
weight loss is due to the loss of two coordinated water-molecules present in the case 
of Fe (HI), Zr (IV) and Th (fV) chelates. On further heating organic component of the 
chelates decompose leading to formation of metal oxide'^. The activation energy was 
calculated from the Broido''* plots of In (In My) vs r' x 10\ Since the plots are 
straight lines, the rate of the reaction is believed to be one over the entire range of 
decomposition. The higher value of activation energy suggests the higher stability. 

The d. c. electrical conductivity of ligand and its polychelates was measured in 
their pellet form in the temperature range 300 - 520 K. (Table 1). The electrical 
conductivity varies exponentially with absolute temperature according to relation a = 
00 exp {-EalkT), where Oo is a constant, Ea is activation energy and k is the Boltzmann 
constant. 

In all the cases electrical conductivity increases with increasing temperature 
indicating these chelates lie in the range of typical semiconductor'^. It was observed 
that, there are two distinct regions. In the low temperature region the slope of plots 
have small values and extrinsic conduction may present in the chelates. In the high 
temperature region, a linear dependence of log o = F (1/7) was observed. In this 
temperature domain the chelates may behave as intrinsic conduction. The lower 
temperature range is the region of extrinsic semiconductors where the conduction of 
carriers from linear localized band to the conduction band. At the higher temperature 
range, the intrinsic region is reached where carriers are thermally activated from 
valence band to the conduction band. The small activation energy observed may be 
attributed to the interaction between the electrons on the cations and n orbitals of 
ligand. This interaction will localize the % electrons charge on the ligand 
molecules'^’'*. 
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Abstract 

The nature of lon-solvent interaction involved in the system of the electrolyte BiBrs in MeCN, 
MeOH and DMF and the effect of addition of EtOH to it at different temperatures using 
electrical conductivity principle, is presented in this paper Investigation showed the 
preferential solvation of the cation by the added ethanol in the case of MeCN+EtOH and 
DMF+EtOH mixtures No preferential solvation by the ethanol is found in the case of 
MeOH+EtOH system The limiting molar conductance of the electrolyte was determined by 
Kraus-Bray and Shedlovsky conductivity models. Limiting molar conductance increased in the 
order of EtOH<DMF<MeOH<MeCN, which is on the expected trend Values of dissociation 
and association constants were also determined The paper also reveals the energy of activation 
involved m the conductance process, thermodynamic parameters of association and solvation 
and the Walden product have been computed From Fuoss-Accascma equation lon-pair 
formation constant Kp and ion triplet constant Kf were calculated and were used to support the 
variation in limiting conductance. Bom relation was also verified An attempt has been made to 
compute the solvation number using mobility method. 

(Keywords : bismuth bromide/solvation behaviour/conductivity/Walden product) 

Introduction 

The measurement of conductance gives important information regarding ion- 
solvation’’^. The study of variation of conductance with concentration, viscosity and 
dielectric constant provides satisfactory information not only about the behaviour of 
solutions but also of ion-solvation^. In recent years, there has been an increasing 
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interest in the behaviour of electrolytes in non-aqueous solvents with a view of 
investigating ion-ion and ion-solvent interactions under varied conditions'*. The 
chemistry of bismuth is diverse^ Inorganic bismuth salts were the first compounds to 
be recognised for therapeutic utility^ Some compounds have been approved for 
human use under Federal Food, Drug and Cosmetic Act for more than 30 years^ 
There are no reports in literature regarding the conductivity behaviour of bismuth 
compounds except for some early investigations^’^. From this laboratory, we have 
communicated the conductance behaviour of Bi(N 03)3 in water + DMF/DMSO*, and 
that of BiBr 3 in MeOH/MeCN + DMF’. As an extension of it, herein, we report the 
solvation behaviour of BiBr 3 in acetonitrile (MeCN), methanol (MeOH), dimethyl 
formamide (DMF), ethanol (EtOH) and their mixtures as a function of temperature. 

Materials and Methods 

Dimethyl formamide (DMF), ethanol (EtOH), methanol (MeOH) and acetonitrile 
(MeCN) were purified as mentioned in the literature''*’". Commercially available 
BiBr 3 (NR Chem make, Mumbai) was used as such without further purification. 
Solution of BiBrs was prepared in pure solvents MeCN, MeOH, DMF and EtOH and 
also in the previously prepared mixtures of known composition (v/v) of 
MeCN/MeOH/DMF+EtOH. Since BiBr 3 undergoes hydrolysis to form BiOBr'^, 
conductivity study is impossible in water. Conductance measurements were made 
with a digital direct reading conductivity meter (model C M. 180, Elico make) and a 
dip type calibrated 8 conductivity cell with cell constant 0.999 cm"'. All the 
measurements were made in a thermostat maintained at the desired temperature 
±0.0 rC. The instrument was standardised as described*’®. 

Results and Discussion 


Limiting molar conductance: 

The solution of BiBr 3 in various compositions of MeCN+EtOH, MeOH+EtOH 
and DMF+EtOH (v/v) were subjected to conductivity measurements at four different 
temperatures (288-318 ± 0.01 K). The conductivity meter gave the specific 
conductance of the solution directly which was used in the determination of molar 
conductance (A„). Debye-Huckel-Onsager (DHO) equation was used to evaluate the 
limiting molar conductance at the initial stage. But the plot of vs Vc was non- 
linear and therefore did not proceed further to analyse the data by DHO equation. 
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Hence, Kraus - Bray equation related to 1:1 electrolyte was used since BiBrs 
undergoes ionisation in a solvent in the following manner'^ 

2BiBr3 = [BiBrz]^ + [BiBr4]- (1) 


The Kraus - Bray equation’"* for such systems may be represented as 
1 A,C 1 


( 2 ) 


where is the limiting molar conductance, C is the concentration in molar and Kc is 
the dissociation constant. The Kraus-Bray plot was drawn for all the cases by taking 
1/Am along the 7-axis and AmC along the X-axis, and the plots were found to be linear. 
From the intercept and slope of the above linear plot (Fig. 1) the limiting molar 
conductance (A„°) and dissociation constant (Kc) were obtained and are shown in 
Table (1 and 2) respectively. The measurements were made at four different 
temperatures and at compositions of (v/v) MeCN/MeOH/DMF+EtOH (0, 20, 40, 60, 
80 and 100% EtOH). The limiting molar conductance obtained are not absolute as it 
does not account for the effect of ionic mobility and activity coefficient on 
conductivity. Therefore, Shedlovsky model'^ was tried, which accounts for the above 
facts, and hence expected to give absolute limiting conductance. The Shedlovsky 
equation may be represented as. 


1 

SA° 

m 


flSK^ 


K 


02 


- + - 


A' 


02 


where 


S = 


2A.^ 3/2 ^ 


( 3 ) 


- 1 .8246 X 1 O’* {CaT I (e Tf^ . 
1 + 50.24 X 1 0* i? {CafVie Tf^ ’ 


log/± = 
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8.20xl0';i°„ 82.5 

(erf '"Ti(€ry/^ 

a = SXJXj ; R = q = e^IlzkT 




Fig 1 Plot of l/A„ (ohm .cm'^ mol) (A„C) (mho cm’’) for BiBrj in MeCN+EtOH system. 

A. 100% MeCN, B . 80% MeCN, C. 60% MeCN, D. 40% MeCN, E. 20% MeCN 
F. 0% MeCN (100% EtOH), x 288K; A1 1298 K; □ 0308 K, 0 3 18 K. 

whm K is the ass^iation constant and/, is the mean ionic activity coefficient of the 
electrolyte. From the intercept and slope of the linear plot of l/5An. vs CA Sf^ the 
hmmng conducmnce and fte association constant (A.) wem dete^ttae^ 

am shown respectively in Tables 1 and 2. < values obtained by both the models am 

almost the same. 
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Table 1- Experimental molar conductance (A°„mho cm^ mol"') at infinite dilution for BiBrj from 
different models in various solvent compositions (v/v) of MeCN/MeOH/DMF+EtOH at- 


r(K) 


0 

20 



40 


60 


80 

100% EtOH 


1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 






MeCN 4- EtOH 






288 

126 

126 

96 

96 

79 

79 

65 

65 

55 

55 

29 

30 

298 

138 

138 

106 

106 

89 

89 

72 

72 

60 

60 

34 

35 

308 

147 

147 

111 

111 

96 

96 

82 

82 

65 

. 65 

42 

42 

318 

163 

163 

122 

122 

100 

100 

89 

89 

74 

74 

49 

49 







MeOH+EtOH 






288 

115 

115 

77 

77 

67 

67 

54 

55 

49 

49 

29 

30 

298 

128 

128 

91 

91 

83 

83 

67 

69 

56 

' 57 

34 

35 

308 

139 

139 

111 

111 

91 

91 

74 

76 

6! 

61 

42 

42 

318 

147 

147 

125 

125 

105 

105 

83 

85 

67 

67 

49 

49 







DMF+EtOH 






288 

45 

46 

43 

43 

42 

42 

41 

41 

37 

37 

29 

30 

298 

51 

53 

50 

50 

49 

49 

47 

47 

40 

40 

34 

35 

308 

57 

58 

56 

56 

54 

55 

53 

53 

47 

47 

42 

42 

318 

67 

67 

63 

63 

60 

61 

59 

59 

53 

53 

49 

49 


1 . Kraus-Bray Model 





2 Shedlovsky Model 
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Table 2- Experimental values of K, and K, for BiBr 3 in various compositions of 
MeCN/MeOH/DMF+EtOH at- 





10'^ 






Ka 



r(K) 

0 

20 

40 

60 

80 

100% 

0 

20 

40 

60 

80 

100% 






MeCN + EtOH 






288 

16.4 

80 

77 

97 

46 

40 

13 3 

24 6 

39 4 

21 8 

72 3 

61 2 

298 

13 7 

83 

7 1 

93 

54 

57 

11.6 

20 1 

38 9 

22 9 

56 1 

36 3 

308 

13 2 

9 1 

66 

92 

77 

50 

102 

14 1 

30 5 

27 9 

27 1 

34 7 

318 

12 9 

9.2 

94 

89 

73 

50 

97 

13.7 

25 8 

30 2 

22 6 

29 7 






MeOH+EtOH 






288 

29 

5.4 

59 

63 

29 

4.0 

76 3 

33 1 

29 6 

60 4 

66 

61 2 

298 

32 

45 

3.9 

42 

32 

57 

75 1 

43 4 

41 3 

65 1 

77 

36 3 

308 

35 

38 

34 

40 

38 

50 

64 3 

61 7 

48.6 

70.9 

62 

34 7 

318 

4.1 

3.4 

37 

4 1 

44 

50 

43 5 

67 6 

61 5 

75 6 

48 

29 7 







DMF+EtOH 






288 

61 

7 1 

46 

2.8 

20 

40 

30 1 

21 0 

31.5 

66 4 

1103 

61.2 

298 

56 

68 

44 

27 

2.4 

5.7 

33.4 

29 4 

34.0 

64.5 

78 9 

36 3 

308 

52 

62 

42 

3 1 

26 

5.0 

34 3 

31.8 

41 1 

61 5 

78 5 

34 7 

318 

4.9 

5.8 

4 1 

3.6 

27 

50 

37 5 

31 2 

36 7 

43.2 

62 4 

29 7 

From Table 1, 

it is 1 

clear that temperature enhances 

the conductivity for all 


compositions of MeCN/MeOH/DMF+EtOH. As expected, increased with 
increase in temperature in all the cases due to increase in thermal energy and mobility 
of ions. Inspection of Table 1 reveals that the limiting molar conductance in pure 
solvents varies in the order MeCN > MeOH > DMF > EtOH. According to Walden'^ 
conductivity and viscosity are inversely related, i.e. higher the conductivity lower will 
be the viscosity. The viscosity of these solvents varies in the reverse order as 
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expected. Size of the solvent molecule is another reason for such variation in 
conductivity. The conductivity data clearly tells that acetonitrile is a better solvent and 
ethanol is the least preferred solvent as far as the dissolution and dissociation of BiBrs 
is concerned. 

Table 1 also shows that addition of EtOH either to MeCN, MeOH or DMF 
decreases electrical conductivity continuously till the system becomes pure ethanol 
medium. The decrease in conductivity is sharp in the case of MeCN+EtOH and 
MeOH + EtOH, whereas the variation is small in the case of DMF+EtOH. Sharp 
decrease in conductivity is observed from 0 to 20% EtOH in the case of 
MeCN/MeOH+EtOH system. This trend is similar to the one observed in the case of 
addition of DMF to MeOH or MeCN^. This decrease hints at the increased ion-solvent 
interaction and also the solvent solvent interaction in the system. The addition of 
EtOH leads to solvent-solvent interaction forming increased solvent mixture 
molecular size and which is involved in solvation of the cation. Hence here we find 
both the ion-solvent .and solvent-solvent interactions. The magnitude of solvent- 
solvent interaction can be determined by its density and viscosity value. They have 
been measured at different temperatures and reported in Table 3. Experimental density 
and viscosity values of pure solvents are in agreement with the literature values’* 
which are shown in bracket in Table 3. As is evident from Table 3 the viscosity 
increased with the increase in % compositions of EtOH and decreased with 
temperature in all the three cases of solvent mixtures. The presence of electrolyte will 
increase the viscosity further for any of the solvent mixtures or that of pure solvents. 
Therefore, conductance decreases from one composition to another at a given 
temperature. Sharp increase in viscosity is observed from 80% to 1 00% EtOH leading 
to sharp decrease in conductivity. From Table 1 it may also be inferred that EtOH 
preferentially solvates the ionic species of the electrolyte. In the case of DMF + EtOH, 
the viscosity decreases with increase in temperature, whereas the conductivity 
increases. But very marginal variation is observed from 0 to 80% EtOH indicating the 
least solvent-solvent interaction. But after 80% EtOH (which is ethanol rich region), 
the conductance decreases indicating higher ion-solvent interaction. That means, here 
it appears either no or negligible amount of solvent-solvent interaction. Probably, 
ethanol does not preferentially solvate the ionic species. The increased ion solvent 
interaction is found only after 60% EtOH. 

Dissociation /Association Constants : 

2BiBr3 = [BiBra]^ + [BiBr4]' 
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Tables- Experimental values of density {d ■ gtn/cm’) and viscosity (n milli poise) for 
MeCN/MeOH/DMF+EtOH at- 


%EtOH 

288K 


298K 

308K 

318K 


d 

T1 

d 

n 

d 


d 






MeCN+EtOH 





0 

0 7940 

4 183 

0 7899 

3 652 

0 7715 

3 398 

0 7607 

3 225 





(3 55) 





20 

0 7880 

4 497 

0 7800 

3 933 

0 7683 

3 707 

0 7600 

3 400 

40 

0 7898 

4 767 

0 7819 

4 107 

0 7705 

3717 

0 7603 

3 422 

60 

0 7961 

6.290 

0 7870 

5.126 

0 7740 

4 627 

0 7645 

4 081 

80 

0 7975 

9.102 

0 7875 

7 461 

0 7771 

6 117 

0 7692 

5 636 

100 

0 8005 

14 200 

0 7934 

11400 

0 7831 

9 446 

0 7739 

7 655 





MeOH+EtOH 





0 

0 7856 

6.244 

0.7865 

5 422 

0 7772 

4 742 

0 7670 

4 174 



(6 25) 


(5 50) 


(4 8) 


(4 15) 

20 

0 8014 

7 652 

0.7916 

6 321 

0 7821 

5414 

0 7725 

4 851 

40 

0 8039 

9 175 

0 7969 

7 535 

0 7868 

6 437 

0 7784 

5 703 

60 

0 8074 

11600 

0 8019 

9.183 

0.7919 

7 890 

0 7835 

6 643 

80 

0 8153 

14 040 

0 8091 

11.300 

0 7981 

8 873 

0 7892 

7 641 

100 

0.8007 

14 200 

0.7934 

11400 

0 7831 

9 446 

01139 

7 655 





DMF+EtOH 





0 

0.9559 

9 440 

0,9443 

7 937 

0.9361 

7.167 

0 9283 

6413 





(7 960) 





20 

0 9298 

9.490 

0 9190 

7.966 

0.9084 

7 241 

0 8998 

6 499 

40 

0.8975 

9.650 

0 8895 

8 037^ 

0.8784 

7.278 

0.8703 

6.559 

60 

0 8665 

9.985 

0 8581 

8.384 

0 8460 

7 365 

0.8372 

6 660 

80 

0 8310 

11.400 

0 8248 

9 099 

0.8149 

7 949 

0 8070 

7 349 

100 

0 8007 

14 200 

0.7934 

11400 

0.7831 

9 449 

0 7739 

7 655 





(11.30) 


(9.30) 
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The electrolyte BiBra dissociates as shown in above equation. The calculated 
dissociation constant {Kc from Kraus-Bray) and association constant {Ka from 
Shedlovsky) values are shown in Table 2. Kc value decreased with the increase in 
temperature from 0-40% EtOH in DMF + EtOH system and increases there onwards. 
At 100% EtOH, it remained almost a constant. In the case of MeCN + EtOH system, 
the decrease in Kc is found at all compositions expect at 20% and 80% EtOH with 
temperature. In MeOH+EtOH system increase is found at 0 and 80% EtOH and 
decrease at other % compositions with increase in temperature. It may be proposed 
that in some cases the system behaves like endothermic and in some other cases like 
exothermic process. Kc decreases from 0 to 80% except at 20% and 100% EtOH in 
DMF+EtOH system. The Ka increases from 0 to 80% in almost all the cases at a given 
temperature and finally decreased at 100% EtOH. Probably, the self association^ of 
the solvent molecule brings about caged structure in the system capturing the 
conducting species of the electrolyte leading to the decrease in conductance. Similar 
trend is observed in the case of MeCN+EtOH system. The trend in variation oi Kg at a 
given temperature in MeOH-i-EtOH system is quite different. Here Ka decreases for 
the initial addition of EtOH to MeOH (up to 40%) with a further increase in Kg till 
80%. Kg values were used to calculate change in enthalpy, free energy and entropy of 
the association process. The Kg increased and Kg decreased in the order of 
MeOH<EtOH<DMF<MeCN. 


Plot of log Am VS log C was drawn on the basis of Fuoss equilibrium^’, to prove 
that the association is predominating over dissociation process, which expects a slope 
of the linear plot around -0.5. In the present case slope value is more than -0.5. 
Therefore Fuoss-Kraus equation'* (eqn. 4) was used to prove further the presence of 
ion-pair association, where g(c) is a factor which incorporates all inter ionic 
interaction terms K^ and Kj are the 


ACg(C)VC=^ + i^^^ 


A 


‘■m2 


ion pair and triple ion pair formation constants and (a°„)^ is the limiting molar 
conductance of the triple ion. The values of ion pair formation constant Kj, and ion 
triplet formation constant Kr were obtained respectively from the intercept and slope 
of the linear plot obtained from Am VC vs [l-Am/A®j„]C. The value of Kj, is found 
higher than Kr at all compositions at 298K in all the solvent mixtures (Table 4). 
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Moreover, a striking feature in this case is an increase in the value of Kp and decrease 
in the value of iiTr from 20% to 80% EtOH. So the decrease in from 0-100% EtOH 
in all the cases can also be explained on the basis of Kp or Kt variation. This ideal 
trend is not observed in the case of pure solvents. 


Table 4- Experimentally determined values of Kp and Kj for BiBr 3 m MeCN/MeOH/DMF+EtOH 
mixture at 25°C. 




Kp 



‘ Kj 


%EtOH 1 

2 

3 

1 

2 

3 

0 

74 40 

126 50 

125.40 

52 65 

17 53 

40 14 

20 

116.16 

111 71 

108 50 

80 78 

29 04 

41 64 

40 

124.54 

120 21 

131 73 

69 71 

24 84 

27 99 

60 

128 18 

123.71 

138.88 

84.36 

21 00 

20.82 

80 

135.72 

128 52 

151 47 

62 05 

19.52 

18 63 

100 

128 51 

128 51 

128.51 

32.54 

32 54 

32 54 

1 

MeCN+EtOH mixture 

2 MeOH+EtOH mixture 

3. 

DMF+EtOH 

mixture 


Thermodynamic Parameters : 


The energy of activation. Eg of the conducting rate process is calculated from the 
Arrhenius equation^^, Ar^ = where A is the frequency factor, i?, the gas 

constant and T, the temperature on absolute scale. From the slope of the plot of log 
Am vs 1/r, Eg is calculated and are shown in Table 5. The variation in Eg with % 
composition of ethanol in anyone of the systems is almost on the expected trend. In 
the case of pure solvents, conductivity increases in the order, EtOH < DMF < MeOH 
< MeCN, where as the energy of activation increases in the reverse order. That means 
ideality is fully followed. the enthalpy of association is calculated from the slope 
of the linear plot of log Kg vs l/Tin all the cases. It is found to be positive from 0 to 
40% EtOH in DMF+EtOH mixture and 0 to 60% EtOH in MeOH+EtOH mixture 
indicating the endothermic behaviour of the system, but in other compositions 
including MeCN+EtOH, the AHg is found to be negative and hence the system 
behaves like exothermic. AG^, the free energy of association is calculated using the 
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equation AGa - -RT In ka at different temperatures (the value at 308K is shown in 
Table 5) and in all the cases is found to be negative indicating both the exothermic 
and endothermic behaviour of the system under different conditions is feasible. ASa 
the entropy of association (value at 308K is shown in Table 5) is found to be positive 
in 0 to 60% EtOH in DMF+EtOH and MeOH +EtOH and in rest of the cases it is 
found to be negative. The values of AHa, AGa and are given in Table 5. 


Table 5 Calculated thermodynamic parameter values for BiBrj under varying compositions (v/v) of 
MeCN/MeOH/DMF+EtOH 


Thermodynamic 

parameters 

0 

20 

40 

60 

80 

100%EtOH 




MeCN+EtOH 




Ea(kJ mor*) 

64 

66 

7 1 

7.4 

75 

11.5 

AHa(kJ moT') 

-10 1 

-15.3 

-11.2 

89 

-29.2 

-18.2 

AS(K-'kJ mor')* 

-0 013 

-0 027 

-0 01 

0.058 

-0 069 

-0 03 

-AG (kJ moT')* 

59 

68 

8.7 

85 

8.4 

9 1 




MeOH+EtOH 




Ea(kJ mol"') 

67 

9,2 

8.9 

91 

10.8 

11 5 

AHa (kJ moT') 

66 

178 

176 

56 

-19.1 

-18.2 

AS(K-'kJ mol-')* 

0 057 

0 092 

0 090 

0.053 

-0.027 

-0 029 

-AG (kJ moT')* 

104 

106 

99 

10.9 

10.6 

9.1 




DMF+EtOH 




Ea(kJ mol"') 

9 1 

9.0 

9.1 

8.9 

10.8 

11.5 

AHa (kJ mol-') 

62 

14.5 

10 1 

-4.8 

-14.5 

-18.2 

AS(K-'U moT')* 

0 05 

0 075 

0 063 

0.019 

-0.011 

-0 029 

-AG (kJ moT')* 

91 

8.9 

9.5 

10.5 

11.2 

9.1 

* at 308K 
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Walden Product : 


Walden product'’ relates molar conductance at infinite dilution(A„“) to viscosity 
(Tjo) of the solvent as 



Ze,F 

SnrT 


(5) 


where r is the Stoke's molecular radius. The product is calculated at all temperatures 
and % compositions for all the three cases. Walden product remains almost a constant 
within the experimental error either with increase in temperature or % compositions in 
DMF+EtOH system indicating the constancy in the Stoke's molecular radius. But in 
the case of MeCN+EtOH system, variation in conductivity has not been compensated 
by the variation in viscosity for a given temperature. Hence, Stoke's molecular radius 
does not remain the same. In MeOH+EtOH mixture the variation of Walden product 
is still different. It decreases with the increase in temperature at all compositions 
except at 20% EtOH. This decrease in the product from one composition to other is 
observed till about 40% EtOH with a later increase at all temperatures except at 3 1 8K.. 
From the equation 5 it is clear that the plot of Walden product vs \/r, should be linear. 
This is true in the case of MeCN/MeOH+EtOH mixture. 


Corrected Stoke's radii was calculated*’’ using equation 6 in all the solvent 
mixtures and the resulted values are given in Table 6. 


r, = 


ZF^ 

6nNkUo 


0.0103 + r, 


( 6 ) 


where rj, - 0.85 A° for dipolar unassociated solvents and = 1.13 A° applied to protic 
and other associated solvents, N is Avagadro number, Z is ionic charge, F is faraday 
constant and is the viscosity. There is a small variation in the radius with change in 
temperature as well as composition especially in the case of MeCN/MeOH+ElOH 
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Table &- Corrected Stoke's radius (A°) of the conducting molecular species of BiBrs under varying 
compositions oi MeCN/MeOH/DMF-rEtOH at- 


T/K 

0 

20 

40 

60 

80 

100%EtOH 




MeCN +EtOH 




288 

2.79 

3 11 

3 36 

3 45 

3 04 

3 30 

298 

2 86 

3 16 

3 42 

3 64 

3 22 

3 40 

308 

2 86 

3 19 

3 47 

3 59 

3 48 

3 46 

318 

2 77 

3 16 

3.56 

3 67 

3.36 

3.56 




MeOH+EtOH 




288 

2 63 

2 87 

2 80 

2 71 

261 

3.30 

298 

2 65 

2 88 

2.75 

2 72 

2 67 

3 40 

308 

2 69 

2.80 

2 82 

2 76 

2.92 

3 46 

318 

2.76 

2.77 

2 77 

281 

2 94 

3.56 




DMF+EtOH 




288 

3.24 

3 20 

3 13 

3.36 

3 35 

3 30 

298 

331 

3 25 

3 26 

3.43 

3 66 

3 40 

308 

331 

3 30 

3.23 

3 50 

3.61 

3 46 

318 

3.24 

3 24 

3 23 

3.49 

3.50 

3.36 


Thermodynamics of Solvation : 

Thermodynamic parameters, enthalpy free energy (AG^-s) and entropy 

of solvation have been calculated using Bom relation* 


AG,. 


r, 


ej 




^Ss 


N{Z,eof 1 5e 
2r, e** Sr 


( 8 ) 
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= AG,., + TA^w. (9) 

where Na is Avagadro Number, Z, is the charge of the ion, eo, the electronic charge 
and r, is the corrected Stake's radius. AH,., calculated is found to be negative (Table 7) 
indicating the exothermic behavior of the solvation process. AG,., is found to be 
negative indicating the feasibility of the solvation process. The decrease in the value 
of AGs., with increase in temperature indicates that the species prefer to remain at 
lower temperature, where the stability of the species is higher. AG,., decreases with 
the increase in EtOH concentration in all the three cases of solvent mixtures in the 
positive direction, i.e. species try to be either in pure solvents or in solvent rich 
region. The plot of -AG,., vs 1/r, was found linear indicating the verification of 
Bom relation. AiS,., is also negative but very small compared to AH,., or AG,.,. It 
decreased with the increase in the amount of ethanol at a given temperature or with 
increase in temperature for a given composition. That means, the disordemess and 
electrostriction increases with the increase in the amount of EtOH either to MeCN, 
MeOHorDMF. 


Solvation Number - So : 

The solvation number S„ of Bi^"^ ion was calculated in four pure solvents at all 
temperatures using the relation, 







( 10 ) 


r, is the Stake's corrected radius, rcry is the crystallographic radius and r,oi is the radius 
of the solvent molecule. Very for Bi^^ was obtained from the literature^®. The solvation 
number was found to be nearly one for all the solvents except in ethanol where 
it is nearing two (Table 8). So, species became bulky due to which mobility 
decreased and hence conductivity decreased from other solvents to ethanol for BiBra 
species at a given temperature. This is what has been observed in the present case 
(Table 1). 



Table 7- Thermodynamics of solvation for BiBrj in various compositions of MeCN/MeOH/DMF+EtOH at 
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Tabic 8- Solvation number (5„) of the cation of the species of BiBrs m MeCN, MeOH, DMF and EtOH 


at- 


Sn 

Solvent 

288K 

298K 

308K 

318K 

MeCN 

0.95 

1 03 

1 03 

0.93 

MeOH 

1.01 

1 04 

1.09 

1.19 

DMF 

1 06 

1 12 

1 12 

1.06 

EtOH 

140 

1 54 

1 62 

1 77 
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Abstract 

A method for the determination of potassium is reported Ferroin is used as a reagent in 6-10 N 
sulphuric acid medium which gives blue colour with potassium salt solution and the blue colour 
obtained has got 402 nm maximum absorbance where potassium is estimated colonmetrically 
in the range of 3 8-13 25 mg respectively 

(Keywords : colorimetric determination/potassium/ferrom) 

Introduction 

Several methods of determination of potassium have been reported mostly 
involving the instrumental methods viz., nephlometric, gasometric, flamephotometric 
and spectrophotometric nature. The widely used flame photometric method requires 
special apparatus, sensitive detectors, photomultipliers and suffers the disadvantage of 
the interference of other elements which is sometimes contradictory in the 
determination of potassium. 

The method for determination of potassium reported in this paper, is a simple and 
direct one. The potassium bromate in solution on reaction with ferroin in 6-10 N 
sulphuric acid medium gives blue colour and the optical density with reference to a 
known standard solution is measured and the amount is calculated. 

Materials and Method 


Reagents: 

About 0.1 N solution of potassium bromate is prepared by weighing the required 
amount of Analar potassium bromate (supplied by B.D.H. Ltd., England) which has 
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been previously dried at 120°C for two hours and making up the solution to a known 
volume. The solution is standardised iodometrically’ . 

0.025 M ferroin is prepared by dissolving 1.485 g of ortho-phenanthroline mono 
hydrate in 100 mL of 0.025 M ferrous sulphate^. 

All other reagents used are of 'AnalaR' grade. 

Colorimetric determination 

Desired volume of 1 : 1 sulphuric acid is taken in each 25 mL volumetric flask to 
give an overall concentration of 6N. Different volumes of potassium bromate are 
added and mixed with 0.05 mL ferroin reagent. The blue colour developed is stable at 
room temperature. The maximum absorbance of the blue colour developed is at 402 
nm. The Beers law curve obtained from known standard nitrate solution is given in 
Fig. 1. Typical results are given in Table 1. 


X= 402 nm 



Fig. 1-VaIidity of Beer's law - Colorometnc determination of potassium with ferroin 
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Table l-Colonmetnc determination of potassium with ferroin Amount of potassium as K 2 O, mg 


Taken ^ 

Found 

3 82 

3 804 

7 361 

7 358 

10 645 

10 648 

13 250 

13 258 


Results and Discussion 

It is found that the acidity should be in between 6-10 N for sulphuric acid for an 
immediate reaction. Below 6 N the colour development is rather slow and above ION 
the colour developed gets slowly discharged. It is observed that the colour attains 
maximum intensity within 3 min after mixing the reagent and remains stable for 48 h, 
if the reagent is in excess. In the colorimetric determination, the blue colour developed 
is stable at room temperature. It is observed that the blue colour developed has got 
maximum absorbance at 402 nm. The Beer's law is obeyed upto 13.2 mg of K 2 O. 
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Abstract 

The paper deals with the combined harvesting of a prey-predator fishery with low predator 
density, the prey species obeying the logistic law of giowth The existence of a nontrivial 
steady state of the dynamical system describing the fishery is discussed Both local and global 
stability of the steady state are studied It is also examined whether the system possesses any 
limit cycle Conditions for the existence of a bionomic equilibrium are established The optimal 
harvest policy is discussed from the view point of control theory The results are illustrated 
with the help of a numerical example. Sensitivity of the optimal solution to changes in different 
parameter values is then studied by using the same numerical example 

(Keywords Harvesting/prey~predator fishery/bionomic equilibnum/optimal harvest policy) 

1. Introduction 

The problems related to harvesting of multispecies fisheries have drawn attention 
of researchers from time to time. It is always difficult, both theoretically and 
practically, to determine an optimal harvesting policy for the exploitation of a 
multispecies fishery. The difficulty actually lies in the construction of a realistic 
model of a multispecies system which leads to an analytically tractable optimal 
control problem. Carrying out dynamic optimization in problem involving more than 
two state variables is mathematically quite difficult. Moreover, reliable estimates of 
interaction coefficients are available for very few multispecies communities. Inspite of 
these difficulties, researchers are engaged in the study of multispecies fisheries from 
mathematical view points. Mercer' made a survey of multispecies approaches to 
fishery management. Clark^ devoted a separate chapter to this topic in his first book. 
In his second book, he discussed the prospects of bionomic modelling in the sphere of 
multispecies fisheries. Acknowledging the difficulties in multispecies modelling. 
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Clark^ observed: “It should be clear that no amount of bionomic modelling is going to 
solve the problems of multispecies fisheries. The best that can be hoped for is that a 
useful classification of the issues, and perhaps some insights into appropriate 
management principles will emerge.” Clark' discussed an optimal equilibrium policy 
for combined harvesting of two ecologically independent species. 
Mesterton-Gibbons'' extended this model to discuss the optimal approach to 
equilibrium. Chaudhuri^'^ discussed a bionomic model of nonselective harvesting of 
two fish species competing for an external resource supporting their growth and 
survival. Saha Ray and Chaudhun’ discussed a Lotka-Volterra prey-predator model 
with harvesting and environmental perturbations. Bionomic exploitation of a 
Lotka-Volterra prey-predator system was studied by Chaudhuri and Saha Ray*. The 
optimal policy for combined harvesting of a prey-predator community was also 
discussed by Mesterton-Giobbons®. In a recent paper, Mesterton-Gibbons'^ 
reconsidered the problems of Chaudhuri*'® to derive its optimal solution. 

In the present paper, we have adopted a model similar to that of Smith" to 
describe the development of a prey-predator community in a fishery. Smeach and 
Rust'^ made a stochastic study of the model of Smith" using the mean-square theory 
of stochastic differential equations. However, they were not concerned with 
harvesting and as a consequence, they had nothing to do with bionomic analysis We 
considered joint harvesting of a prey-predator community whose biological 
development is governed by a Smith type model" The dynamical behaviour of the 
system is first studied. It is shown that a nontrivial steady state exists under certain 
conditions which are amenable to interpretations relevant to the model. Both local and 
global stability of this steady state are then discussed; existence of limit cycles is also 
examined. Bionomic equilibrium of the system is next discussed. We then study the 
optimal harvesting policy from a control-theoretical view point. A numerical example 
is taken up to illustrate the results. Sensitivity of the optimal solution to changes in the 
values of diflferent parameters is also examined. 

2. Statement of the Problem 

Let us consider two competing fish population whose growth obey the following 
dynamical system : 


i&i _ 

It 




-oa,JC2 


dx2 _ 

~dt T" 


^Px2 =G(j:,,X2) 


( 1 ) 
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where r\, ri, a, P, k are all positive constants. 

Here Xi(/) = the size or biomass of the prey population at time t, 

X 2 {t) = the size or biomass of predator population at time t, 

r\ = biotic potential of the prey, 

r 2 = the maximum reproductive rate of the predator, 

a = the rate at which the prey species is removed due to predation effects, 

P = per capita death rate of the predator, 
k = the canying capacity of the prey population. 

The crowding effect term (due to intra-specific competition) -5-1- is present in 

k 

the first equation and it is absent in the second equation of the model. This indicates 
that the predator density is low; the predator species dies out exponentially in the 
absence of the prey. 

Equations of the type (1) have been used by Smith" to model the development of 
a prey-predator community with low predator density. 

We now assume that both the species are subjected to a combined harvesting 
effort governed by the differential equations: 




1-i 


- a JCjXj - q\Ex^ 


c&c, 

dt k 2 </2 2 


( 2 ) 


where qy, are the catchability coefficients of the two species and E denotes the 
effort devoted to their combined harvesting. We intend to study the characteristics of 
the exploited prey-predator fish species described by the dynamical system (2) and 
also some economic implications of the harvesting policy. 
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3. Dynamical Behaviour 

fAJ Steady state : 

The steady states of the system of equations (2) are 

/’o(0, 0), P< (j:,*. xl] P2{x„ X,) 


where x*, = — (r, - QiE), 

(3) 

n 

xl = Q, 

(4) 

_ki^ + q2E) 

(5) 

^2 

and X 2 = {ri(r 2 -^)-E{riq 2 + r 2 qi)}/ar 2 . 

(6) 


The nontrivial steady states of the system (2) are 
Xj) andPzfe, Xj). 

For jc,’ > 0 , we must have from (3), that E < — . (7) 

h 

From (5), we have x, is always positive. 

It is obvious that r2>(3 (8) 

and therefore. X 2 > 0 iff ri(r 2 - P) > E(ri ^2 + fjiji) ' (9) 


For the unexploited system (E = 0), the nontrivial steady states P \ and P 2 


become P i(k, 0) and P 2 


0^2 




, f2>P. 
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These results perfectly agree with those of Smeach and Rust’^ for (3 = 1. 
Comparing Pi with P i, we find that the combined harvesting effort incre^es the 
equilibrium level of the prey, but decreases that for the predator. Comparing P \ with 
/*,, it is observed that, in the absence of the predator, the harvesting effort decreases 
the prey level. 

The condition (9) may be written as : 


^ )"■} 

0 < — + — < 


n(^2-P) 

Eqxqi 


( 10 ) 


In fishery literature, the ratio of biotic potential (r) to the catchability coefficient 
(q) of a species is called its BTP (biotechnical productivity)^. Thus, ^ 

prey-BTP and the predator-BTP respectively. Condition (10) thus establishes 

^ as the maximum value of the sum of these two BTP’s. When the sum of 

Eqiq2 

these BTP’s attains the maximum value, the predator population is extinct (x 2 = 0). 
On the other hand condition (7) shows that, in the absence of predators, if the 
harvesting effort is equal to the prey-BTP, the prey population is extinct ( x* = 0) and 
Pi coincides with Pq. 

Condition (9) may also be rewritten as. 


o<£:< 


'1^2 +'*291 


(H) 


Condition (11) defines P„,ax ^ as the upper limit of the level of 

m 2 + ^ 2 ?! 

harvesting effort. We must have 0<E< jBmax for the existence of positive steady states 

f ~~ 13 ') 

of both the prey and the predator. If P = Pmax ~ — I— ^ the predator population 

m2+r2qi 


becomes extinct (x 2 = 0). 
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(B) Local stability 

The variational matrix of the system (2) is 
2r^x^ 


V{xi,x^) = 


n -- 


--00:2 -( 1 \E 


■CXXi 


AX 


2-^2 


^2^1 


•P-92^ 


At Po(0, 0), V(0, 0) = 


r^-q^E 

0 


0 

-i^ + q^EX ■ 


Therefore, the eigenvalues of the variational matrix of the system (2) at Po(0, 0) 
are ri - q\E and -(p + qiE). 

If the prey-BTP — > £ , then the eigenvalue (ri-^i E) is positive while the other 
?! 

eigenvalue is obviously negative. In this case, the equilibrium point Po(0, 0) is a 
saddle point. 

f* 

On the other hand, if — < £ , then both the eigenvalues being real and negative, 

£0 is a stable node. However, — < £ , makes both x\ and 3c 2 negative. In that case, 

no nontrivial steady state exists. We, therefore, discard the case — < £ and accept (7) 

9i 

as the basic condition for the existence of a nontrivial steady state. 


At 


Pi(x,*,x*)i v(xi*,x*)= 


k 


0 



“OCXj 

k J 

(rj-q-i£) 

'l 

— ('i-9i^)-p-?2^ 
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The eigenvalues at P\( x* , xl ) are 

-{rx-qiE) and ^{r^- q^E)-^- q^E . 
r. 


Now -(ri - qiE^ < 0 by (7). 


''l n 

= - h ('-2 - P) - Eiq^r^ + r^q ^ )} < 0, 


then P\ is a stable node. However, ;c 2 becoming negative in this case, P 2 does- not 
exist. 

Thus, P\ is the only nontrivial steady state and it is a stable node when the reverse 
of the inequality (9) holds. When P 2 also exists, the inequality (9) holds and the 
second eigen value becomes positive; as a result Pi becomes a saddle point. 

The variational matrix of the system (2) at P 2 (x 1 , x 2 ) is 


V{xx,x-2) 


The eigen values X,(j = 1, 2) of the system (2) at PaCx i, x 2 ) are the roots of the 
following quadratic equation in X : 


^2X2 
k 

or, kX^ + ri 3c lA, + ar 2 X 1^2 = 0 . 


-OUj 

= 0 . 

-X 



( 12 ) 
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Discriminant of the equation (12) is 
D = -4a . 

V X 

In equation (12), sum of the roots = — ^ < 0 
oc r jc jc 

and product of the roots = — ^ > 0 . 

K 

Therefore, the roots are either both real and negative or complex conjugates with 
negative real parts. 

The eigen values are real or complex depending on Z) > 0 or < 0 which implies 


E>or< 


4^2(^1^2+92n)+?2'l^ 


= E (say) 


(13) 


This result perfectly agrees with that of Smeach and Rust'^ for p = 1 in the case of 
an unexploited system (E = 0). Combining (13) with the condition on for the 
existence of the steady state (3c 1 0, 3c 2 0), we find that 

(a) E ^E < Efnax 

when the eigenvalues are real and negative and 

(b) 0<£<Min(£„ax. £) 

when the eigenvalues are complex conjugates with negative real parts. 

When the eigenvalues of the system (2) are both real and negative, the non-trivial 
critical point ? 2 ( ^ i, ^ 2 ) will be an asymptotically stable node. On the other hand, the 
critical point P 2 will be a stable focus when the eigenvalues are complex conjugates 
with negative real parts. 

Thus, when both Pi and Pj exist, i.e. when both the inequalities (7) and (9) hold 
simultaneously, we find that 

(i) Pa is a saddle point; 

(ii) P\ is a saddle point; 

(iii) P 2 is a stable node or a stable focus depending on the value of the effort E. 
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By this eigenvalue analysis, we get the behaviour of the dynamical system (2) 
only in the neighbourhood of the critical points. 
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fC) Limit cycle : 

The system may have a limit cycle, we now examine the possibility of existence 
of a limit cycle of the non-linear system (2) by using Bendixon-Dulac test in the form 
suggested by Clark^ (p.208). 

Writing the dynamical system (2) in the form 


A, 

dt 

dx2 

Ht 


= ^1 («! + C, ) = /(x, , Xj ) 

= ^2 («2^1 + ^2^2 + ^2 ) = g(x, , X2 ) 


(14) 


where a\ = 



At 


p ^2 = 0 , C2 = - (P + ^2-^). 


We find A = a\b 2 - aibx - — 2.:^0 

k 

and a,c2(^)i - 62) + i2Ci(a2 - 00=-— 

k 

A ’ A r2 ■ 

Introducing the function 

-0-1 

5(X!, X 2 ) = x,‘’"'x2"‘ ^xf’xj'^ , 
wefindthat (p+g2^)<0 

C«X| CfX^^ 
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always in a simply connected region D in the x-x^ plane. 

Hence the expression 

d ^ 

— {Bf) + - — (Sg) does not change sign in D. 
dxi 0X2 

Hence the system ( 2 ) does not possess limit cycles. 

(D) Global stability : 

We now examine the global stability of the dynamical system ( 2 ). We have 
already noted that this system has a unique positive equilibrium at P2{x 3 c 2) 
provided the conditions ( 8 ) and ( 9 ) hold. The prey-density having a carrying capacity 
and the predator density being low, we may assume that the solutions of the dynamical 
system ( 2 ) are eventually uniformly bounded in a finite region B in the (x^, x^) phase 

plane and hence all solutions of ( 2 ) starting in R2 enter into the region B. We shall 
now prove that the interior equilibrium P2( x u 3 c 2) is globally asymptotically stable in 
the Lyapunov sense in the region B. 

We define a Lyapunov function 




r 

(x "ll 

F(jc, ,X2 ) = - X, - Xj In 

Xj 

+ d\x2 -X2-X2 In 

^2 1 
_ 

Jlj 

1 



where c/ is a suitable constant to be determined in the subsequent steps. 

Obviously Fis a positive definite function in the region B except at P2(x i, x 2) where 
it vanishes. 

Further , Lini ^ F(xi, X2) = . Lim V(xi, X2) = «> 


The derivative of V along the solution of ( 2 ) is 
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Clearly V =0 when xi = X]. On the set {(xi, xz): Xi= x 1} the only invariant set is 
{(xijXj)}. Hence by Lasalle’s invariance principle’^ P2(^]>^2) is globally 
asymptotically stable. As a consequence of the global asymptotic stability of the 
equilibrium point P2 in the subspace 5 of R2 , it again follows that the system ( 2 ) has 
no limit cycles. 


4. Bionomic Equilibrium 


Let, c = fishing cost per unit effort, 

= price per unit bio-mass of the prey species 

and = price per unit bio-mass of the predator species. 
The economic rent (net revenue) at any time is given by 
n(xi, X2, E) =piqiExi + p2q2Ex2 - cE. 
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We find that, =0 = X2 yield the results 

(15) 

(16) 

Eliminating E, we find that, the nontrivial equilibrium point of the exploited system 
(2) occurs at a point on the line 

r2Xi arjX, 6 riCo-p) 

*92 fi92+''29l 92 n92+''29l 

The biological equilibrium line (17) meets the a:i - axis and jC 2 -axis at (3q,0) and 
( 0 ,X 2 ) respectively where 

_ *92 f P , ^i('2 -P) ] _ *(p9i +^ 92 ) . Q 

‘ ^2 \q2 fi92+''29lj fi92+^29l 

and - (fi92 +^29i) [ P , fi (^2 "P) 1 _ P9i +'192 ^ q 

^ Ct^2 192 ^l92+f29lj “92 

The bionomic equilibrium (xi«,, X 2 «) of the model (2) is given by (17) together with the 
condition 


^ 2^1 


P 


+ E, 


and 


*92 92 

aj ^2 _ '1(^2 


■P) 


'■i92+^29i n92+''29l 


= (P,9,Xi+P292^2“^)^" ® (18) 

As long as (p^q^x^ + ^ 2 ^ 2^2 ~ 0 points on the equilibrium line (17), the 

species remain unexploited because the fishery fails to produce any positive economic 
revenue. 

The three situations that may arise with respect to the intersections of the lines 
(17) and (18) are given below in corresponding figures. 
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In Fig. 2(a), positive bionomic equilibrium exists for both the species. 

In Figs. 2(b) and 2(d), the zero-profit line and the bionomic equilibrium line 
intersect in the second quadrant of the phase plane where is negative. Hence no 
bionomic equilibrium exists in these cases. Similarly, no bionomic equilibrium exists 
in Figs. 2(c) and 2(e) also in which the two lines intersects in the fourth quadrant 
where X 200 is negative. 

In Figs. 2(b) and 2(c), the biological equilibrium line lies above the economic 
indifference line. Thus there is an economic advantage in both the cases whatever the 
biological equilibrium point is. It is a case of economic equilibrium only which is 
attained before the biological steady state is reached. 

In Figs. 2(d) and 2(e), the economic indifference line lies above the biological 
equilibrium line. As a result, the fishermen would continue to harvest until the 
zero-profit lines rt = 0 is reached without caring for biological equilibrium. Thus, this 
is also a case of economic equilibrium only occurring beyond the biological steady 
state. 

Solving (1 7) and (18), we get. 


A:{ac-p2(r,q-2+P?i)} 

.{a^i5,-P2(fi92+^29i)} 

(19) 

_ kpiqiimi +?iP)-c(^i?2 +''2«i) 
a^,g,?2-P292(n?2+^2^l) 

(20) 


Thus jciM and X 2 a are respectively given by (5) and (6) in which E is to be replaced by 
the bionomic effort level E„. Using the open-access condition %(£„) = 0, we get 

^ _ g cr; - aPpi^ifc - (r^ - p) 

ap^q^q^k-p^giimi+r^q^) 


Now xia, > 0 provided either of the following conditions holds 


(A) p > max 


ac 


^ ] 

/i?2 + P?i’ '■i92+^2?J’ 
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Xi Xi Xioc 

!d) (e) 


Fig. 2- Bionomic equilibrium of the prey-predator fishery 
(a) Non-extinction,(b) and (d) Extinction of prey, (c) and (e) Extinction of predator 
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(B) p < min 


ac 





Similarly, > 0 provided that 


?1 92 1 I 92 j P2^2 J 


(D) —+—> mini- 
mi 92 ^ 




H 

92 J 


|,S^|h0lds. 

^292 , 


When either (A) or (5) and (C) or (D) hold simultaneously, both the species co-exist 
in bionomic equilibrium. When either (A) or (B) holds and neither (C) nor (D) holds, 
the predator species faces extinction in the bionomic equilibrium; the bionomic 
equilibrium of the fishery is supported by the prey species only. Similarly when either 
(C) or (Z)) holds and neither (A) nor (B) holds, the prey species faces extinction in the 
bionomic equilibrium and as a consequence, the predator population also is eventually 
driven to extinction. 


5. Optimal Harvest Policy ; 

Our objective is to maximize the present value J of a continuous-time stream of 
revenues given by 


J = [piq^Xy + Pzq2h. ~ dt (21) 

0 

where 6 denotes the instantaneous annual rate of discount. 

We intend to maximize (21) subject to the state equations (2). 

The Hamiltonian for the problem is given by 
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H = e'®' + P2<l2Xi - + K ] r^x, - ^ - cocj x, - q^Ex^ i 


^X^}W2_ 




( 22 ) 


where X|(^) and A-aCO ^re the adjoint variables. 
The adjoint equations are 

dX, dH 


dt 


'—~~-e-“p,q,E-X,[F^ -q,E)-X,G, 


(23) 


(it UX'y 


(24) 


We consider the optimal equilibrium solution of the given problem so that, we may 
write 




q\Xl 92X2 

Equations (23) and (24) can be written as 


^1 _ 

dt 


= -e' 


2/-,X, r- I ^2 ^2 ^2 

r, \^-ax^-q,E — 


V 


(25) 


(26) 


and 


dk 


dt 




(27) 


Equation (25) implies -■X^-ax 2 -q^E = 0 


(28) 




and 


(29) 
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Using (28) and (29) in (26) and (27) we get 


dk^ __ -5r 






and ^ = -e'®'p292^ + ^i“^i- 
dt 


Differentiating (30) with respect to t and using (3 1 ), we get, 


d Xj dk^ ^ r-^]0.x^X2 _ £^-st | 

dt^ k dt k I '' k 


Auxilliary equation of (32) is 


2 r,x, r.ax.x-, . 

k k 


If m\ and mi be the roots of the quadratic equation (33), then. 


r^x, ^ r2ax,X2 „ 


ffji + /«2 = > 0, m,W2 = — — 

k k 


Therefore, the roots of the equation (33) are both real and positive or complex 
conjugates with positive real parts. 

Case (i) : When both the roots are real and unequal, the general solution of (32) is 

(r) = A^e'^' + + ^e'^' 

' ' 1 


where M, = £ 8 ^ 

I A j 

AT -52 , I h^x^Xi 

' k k 
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and Ai, B\ are arbitrary constants. 

Thus the shadow price remains bounded as t ^ if y4i = 5i = 0 and then 

(t) = — = constant. 

Case (ii) : When the roots of the equation (33) are real but equal, each of the roots is 
2k 

The general solution of (32) is therefore, 

where A 2 and B 2 are arbitrary constants. 

In this case also, the shadow price ^ie*‘ remains bounded as t -> 00 if ^2 = ^2 = 0 

and e^'X, (t) = — = constant. 

^ iVi 

Case (hi) : When mi and m 2 are complex conjugates with positive real parts, let mi = 
tti + /Pi, m 2 = ai - iPi where ai > 0. 

The solution of (32) then becomes 

X,,(r)e®' = cos(Pit + s)+ ^ where A 2 and s are arbitrary constants. 

As t — > 00 the first term becomes unbounded in an oscilatory manner. Hence the 
shadow price e®‘ ^i(t) remains bounded as t 00 iff ^2 0 so that we have 

e^’XAt) = = constant. 

iVi 

In a similar manner, we can show that, X 2 (t)e^‘ constant 
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„ r,x,p-,gjE 

where AE = 5 p^q^E - a x,p ^qE + ^ . 

liras the shadow prices e^‘ ?,,(/) ( 2 = 1 , 2 ) remain constant over time in optima! 
equilibrium when the}' satisfy the transversibility condition at infinity'"*’ i.e , they 
remain bounded as / -> co. Now the last condition that the Hamiltonian H must be a 
maximum for an admissible control E leads to the condition 


BE 




+>-2?2^2" 


d% 


(34) 


This shows that the total user cost of harvest per unit effort must equal the discounted 
value of the future price at the optimal equilibrium level. 

Equation (34) can also be written as 


or, 


M] ) M, 


f 




5/2,9, + 


V 


iV, 


ViPi<li 


P\ 




= c 


+ 9,^2 


p ^ f£ E 

2 riX,5 r2X^X2a 


(35) 


Eliminating E from (28) and (29), we get. 


?2 


r, -■ 


r,x 


im 


■x-,a 


= ?i 


-P • 


(36) 


Equation (35) and (36) together determine the optimal equilibrium populations = x 


^ndx^ = x^. 


15 



harvesting of a PREY-PREDATOR FISHERY WITH LOW PREDATOR DENSITY 461 


6. Numerical Example 

Let us take the parameters of the system as 

r^ = Q.l\,r^= 1.8, A: = 600,a = 00009, p = l,c= 100,P| = 250,j?2 = 500, 

8 = 0.15, = 0.0006, q^ = 0.0007, £■ = 20 in appropriate units. 


For the given paramedic values, we find that the biological equilibrium occurs at the 
point (338, 331) and E « 306.26, 360.18. 

Thus 0 < £■ ( = 20)< min ( E , ^max)- 

Therefore, criterion (b) determines that the non-trivial equilibrium point is a stable 
focus. 

For the parametric values prescribed above, it is seen that the criteria (A) and (C) are 
satisfied implying that > 0 and > 0. Computations using (19) and (20) actually 

yield « 389 and « 1 19. It is interesting to note that the bionomic equilibrium 

for the prey species is established at a level higher than its biological equilibrium. On 
the other hand, the bionomic equilibrium for the predator occurs at a level which is 
much lower than its biological equilibrium. The optimal equilibrium solutions for the 
same numerical example obtained by solving (35) and (36) are xjg = 385 and X 2 s = 
134. Thus the optimal equilibrium solution lies between the biological equilibrium 
and bionomic equilibrium. 

For different initial values of jri(/). Fig. 3 depicts the time-course display for the prey 
density :)i:i(t). The population density A:i(t) is seen to undergo an oscillatory behaviour 
for smaller values to time t, whatever be the initial values .rj(O). As t gets large, the 

population stabilizes at the biological equilibrium level. 

Fig. 4 shows the lime-course display for the predator density x,(0 for different initial 

values ;!:2(0). Here also X 2 {f) exhibits an oscillatory behaviour for small times. 
Ultimately, it stabilizes to its biological equilibrium level. 

The phase-plane trajectories of the dynamical system (2) are shown in Fig. 5. Clearly 
the biological equilibrium point P is a stable focus. 




harvesting of a prey-predator fishery with low predator density 463 



Fig 5-Phase plane trajectories of the system represented by equations (2) 

7. Sensitivity Analysis 


Using the above numerical example, the sensitivity analysis is performed to study 
the effects of changes of parameter values on the optimal solution. 

Table I- Effects of changes m the parameter values on optimal prey population and predator population 


Parameters 

% change in parameters 

% change in 

% change in x 

25 


+30 

5 33 

-6 87 

n 

+20 

3 67 

-4 72 


-20 

-4 22 

5 38 


o 

1 

-6 58 

8 36 


+30 

-18.67 

23 63 

^2 

+20 

-13 28 

16 83 


-20 

18 11 

-23 25 


-30 

29.90 

-38 60 


Table! Contd, 
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Table 2 Contd 



+30 

34 82 

-44 59 

k 

+20 

22 92 

-29 27 


-20 

-21 84 

27 71 


-30 

- 32.39 

41 05 


+30 

-2 67 

1 60 

a 

+20 

-1 76 

0.98 


-20 

1 70 

-0 57 


-30 

2 53 

-0 65 


+30 

- 6.30 

69 49 

c 

+20 

-4 06 

48 88 


-20 

4 15 

-49 95 


-30 

6 32 

-75 98 


+30 

3 71 

-44 59 

p > 

+20 

2 43 

-29 27 


-20 

-2 30 

27 71 


-30 

-3 41 

41 06 


+30 

1 96 

-23 55 

Pz 

+20 

142 

-17 09 


-20 

-2 21 

26 60 


-30 

-3 86 

46.42 


+30 

0 08 

-0 97 

6 

+20 

0 05 

- 0.66 


-20 

-0 06 

0 70 


-30 

-0 09 

1.06 


+30 

0.23 

-36 70 


+20 

016 

-24 36 


-20 

- 0.18 

24 09 


-30 

- 0.28 

36.05 


Table 2 Contd 
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Table 2 Contd 



+30 

5 08 

-25 83 

^2 

+20 

3 50 

-18 76 


-20 

-4 02 

28 88 


-30 

-6 27 

49 88 


+30 

-0 28 

341 

E 

+20 

-0 19 

2 27 


-20 

0 19 

-2 24 


-30 

0 28 

-3 35 


+30 

19 88 

-25 40 

p 

+20 

13 25 

-16 90 


-20 

-13 25 

16 82 


-30 

-19 87 

25 21 


From the sensitivity analysis in Table - 1, the following points are observed. 

(1) both and are slightly sensitive to changes in the parameter r^, the 
biotic-potential of the prey. increases while decreases with increase 


(2) Both Xjg and x^^ are moderately sensitive and x^^ is more sensitive than jJCjg to 
changes in r 2 , the biotic-potential of the predator. Also decreases while 
x^, increases with the increase of /"a. 

2o 

(3) Both Xjg and x^^ are quite sensitive and x^^ is more sensitive than to 
changes in k, carrying capacity of the prey. As k increases, increases but 
x^, decreases. 

Zo 

(4) Both and are insensitive to changes in a, the predation coefficient. 

10 25 

(5) Xjg has very low sensitivity while possess very high sensitivity to changes 
in c, the cost per unit effort. As c increases, x^g decreases slowly, but X 25 
increases rapidly. 



466 


K. S. CHAUDHURI and T. PRADHAN 


( 6 ) jCjg is slightly sensitive while is highly sensitive to changes in pi, the 
price per unit biomass of the prey. As p] increases, Xjg increases slowly, but 
x^^ decreases rapidly. 

( 7 ) Xjg is almost insensitive while x^^ is quite sensitive to changes inp 2 , the price 
per unit biomass of the predator. :ic,g increases slightly but x^^ decreases 
considerably as p 2 increases. 

(8) X and ;c are almost insensitive to changes in 5, the discount rate. 

( 9 ) JCjg is almost insensitive while x^^ is considerably sensitive to changes in qi, 
the catchability coefficient of the prey. As qi increases x^^ decreases and 
vice-versa. 

( 10 ) has low sensitivity and x^g has high sensitivity to changes in q 2 , the 
catchability coefficient of the predator, undergoes slow increase and x^^ 
undergoes rapid decrease as q 2 increases. 

( 11 ) jCjj is almost insensitive and x^g is slightly sensitive to changes in E, the level 

of harvesting effort. It should, however, be remembered that, this result 
holds for admissible values of ^jas per conditions (a) and (b). 

(12) Both jCjj and x^^ are quite sensitive to changes in (3, the percapita death rate 
of the predator, increases and x^^ decreases with the increase of p. 


8. Discussion 

The present model is suitable for a fishery having a prey-predator community in 
which the prey density is high and that of the predators is low. An appropriate 
example might be a krill-whale fishery. For the numerical example taken here, we 
have determined three equilibria; biological equilibrium jPi( 338,331), bionomic 
equilibrium P2(389,119) and optimal equilibrium P3(385,134). Pi is the most 
conservative of these three solutions. However, the economists may dislike it because 
it takes no economic factors into consideration. Although P 2 is based on economic 
considerations, its acceptability also is questionable because the fishery is, in this 
case, run on the no-profit-no-loss basis. P 3 seems to be the best choice for the fishery 
manager. 
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As far as the dynamical behaviour of the system is concerned, the existence of 
positive steady states for both the species depends primarily on their BTP’s . These 
BTP’s, on the other hand, depends on what types of fishing gear are being used, 
technical efficiency of the fishing gear is reflected in the catchability coefficients q\ 
and q 2 - Thus the choice of the fishing gear is crucial for the survival of the species. 
The BTP’s in turn, define upper and lower limits of the level of fishing effort E. These 
limits of E determine the behaviour of the steady state. Non-existence of any limit 
cycle of the system implies that the system does not undergo periodic oscillations. As 
we have seen in the case of local stability of the system, there are strict eigen -value 
inequalities : either (i) X] < 0, A .2 < 0 or (») Xi, Xi are complex and Re A,, < 0 (/ = 1, 2). 
These strict eigenvalue inequalities are preserved by small changes in the elements of 
the eigen matrix. Thus the stable node in case (/) and the stable focus in case (//) are 
stable relative to small perturbations of the parameters ri, a, p and k. Hence the 
dynamical system (2) represents a reasonable model of a biological system. The 
global asymptotic stability of the steady state implies that the behaviour of the 
solution trajectories does not change drastically if the model is subject to small 
perturbations. 

As we have already seen in Figs. 2(b) and 2(d), practically no bionomic 
equilibrium exists because (17) and (18) intersect at a point with xi«, < 0. In such a 
case, as Clark^ (p. 305) says, the two lines do not intersect, and bionomic equilibrium 
occurs at the poipt (0, X 200 ). This view is untenable because the bionomic equilibrium 
state is defined as the point of the state plane in which both the conditions (17) and 

(18) are satisfied. It is also seen from (19) that xi„ 0 and xin -- ---- ( a finite 

q^a. 

quantity) only when p 2 is finite and p\ -> 00 . This cannot be a practical situation. The 
same observation holds for the situations depicted in Figs. 2(c) and 2(e) also. Thus 
there cannot be any bionomic equilibrium with one species extinct as concluded by 
Clark^. Positive bionomic equilibria of both the species exist when p\ and pi are both 
finite and other conditions, as discussed in section 4, are satisfied. 
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Abstract 

Creep transition in a circular non-homogeneous cylinder, subjected to torsion has been studied 
by transition approach Non-homogeneity in the material is considered due to radial variation 
of modulus of rigidity. Stresses for steady state creep have been derived and discussed 
numerically Effects of the non-homogeneity appear to be prominent and the same shear stress 
distribution has been found for different combinations of the non-homogeneity parameter k and 
the measure n Moreover, non-homogeneity of a certain type having more value at the outer 
radius than at the centre of the cylinder may be beneficial from a design point of view 

(Keywords : creep/torsion/generalised measure) 

Introduction 

The 327 page paper presented to the French Academy by St. Venant\ about one 
and half century ago remains the most outstanding work on the torsion of cylinders of 
arbitrary cross-section, that, is circular, elliptical, rectangular, equilateral, triangular 
etc., and the topic entered many texts^'^ Lee^ considered the torsion of an elastic 
perfectly plastic cylindrical; bar and showed that the curvature of the shearing stress 
lines in the plastic zone has an upper bound proportional to the twisting angle, while 
Cho and Findley^ studied creep and plastic strains under side steps to tension and 
torsion for 304 stainless steel at 593°C. An analysis for the creep deformation of a 
thin-walled cylinder subjected to axial force, bending and twisting moment, was 
carried out by Lau et. al} under the assumptions of Prandtl-Nadai creep law and 
Incompressibility of the material. 

The problem of torsion of a semi-infinite isotropic, homogeneous, elastic solid 
when a circular cylinder welded to its plane boundary and is forced to rotate about its 
axis was first discussed by Reissner and Sagoci^. Erguven'® worked on the Reissner- 
Sagoci problem for a transversely isotropic, non-homogeneous elastic solid by varying 
the modulus of rigidity as a power of the radial co-ordinate and a note on the twisting 
of a non-homogeneous sphere has been furnished by Chaudhuri and Datta^^ on the 
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basis of classical theory of elasticity. In an investigation Cowin'^ observed that for a 
certain class of elastic cylinders with shape intrinsic orthotropy, the solution to the 
torsion problem is the same as the solution to the torsion problem for the isotropic 
cylinder of the same shape if the isotropic shear modulus were replaced by the 
orthotropic cylinders by considering the finite strains incompressibility of the material 
and Norton’s law of creep. 

In this paper an attempt has been made to study creep transition in torsion of a non 
homogeneous cylinder with the help of transition approach. Transition theory 
developed by Seth''*’'^ does not require any adhoc assumptions like incompressibility 
of the material or creep strain law, and solves a more general problem from which 
cases pertaining to such assumptions can be worked out'®'^^. It utilises the concept of 
generalised strain measure and asymptotic transition through the critical points of 
the differential system defining the deformed field. For uni-axial case, Seth'"* has 
defined the generalised strain measure as 

c = ( l/«{ l-{/o//}T (1) 

where n is the measure, m the irreversibility index, 4 and / are the initial and strained 
length of the rod respectively. 

Here, a secondary stage of creep, i.e., m= 1 has been considered and non- 
homogeneity in the material is taken due to radial variation of modulus of rigidity 

H = II or'* (2) 

where po and k are real constants. 

Governing Equations 

Considering a long circular cylinder of radius a subjected to finite twist. Let r, 0, z 
refer to the radial, trangential and axial directions of the cylinder and z-axis is 
coincident with its long axis. The components of displacement are given by 


M= r ( 1 -(3), v= a r z, w = d^z, (3) 

where d\ is a constant, a is the angle of twist per unit length and P is a function of 
only. 
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The generalised components of strain from eqn. (1) are 

err= l/« +P)"], eee = l/« [ 1-P" ], 

= 1/n [c/ 2 - (a^P)"], ear = 1/n (ar)"'^ p”, !► (4) 

^rB ~ ^zr ~ Oj 

where 

c /2 = 1- (1-c/i)" and p' = c/p/c//- 
Stress-strain relations are 

Ty = + 2p,ey, (/,y= 1,2,3) (5) 

where Ty and gy are the stress and strain tensors respectively, "K and p are Lame’s 
constants and 5y is Kronecker’s delta. 

Substituting eqn. (4) into eq. (5), one gets 

Xrr = "kSx + 2p [ 1- ( rp' +P)"] / n, tee = 'kSx + 2p [ 1-p"] / n, 
t„= XSx + 2p [ c /2 - ( arp)" ] / «, tez= 2pp" ( ar)”^V n, T (6) 
t/*e tez 0, J 

where S\ = [2+c/2-P"-(a'P'+ P)'’-( arp)"]/«. 

Equation of equilibrium to be satisfied is 

Strr/5r+(trr -tee)/)^=0 (7) 

Substituting from eqn. (6) into eqn. (7), one gets a non- linear differential 
eqn. in P as 
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(X+ 2 n)pP(P+l r'rfP/^/p=2n[ 1-(P+ 1 )"]/«+2fa'r[p'"-(P+ 1 )7« 
-X{ary-P[X (l+aV)+ (>.+2^) (P+1)”], (8) 

where rp =pP and p. = d]ildr 

Transition points of P in eq. (8) are 


P-> 0, -1, ± 00. 


Torque applied is 


XQzdr (9) 

0 

Statement and Mathematical Model of the Problem 

Following the transition approach, we analyse the problem of creep transition in 
torsion of a cylinder made of the material having variable modulus of rigidity along 
the radius, and governed by eqn. (2), (6) to (9). It is of interest to investigate the 
actual variation of the creep stresses on the assumption that the radial stress vanishes 
on the free surface and no longitudinal force acts on transverse sections of the 
cylinders i.e. 

Trr = 0, at r = a (10) 

a 

and |2nr x^^dr = 0 (11) 

0 


Solution of the Problem 

It has been shown by several workers that the asymptotic transition through 
the principal stress difference gives creep stresses at the transition point P-¥-\. 
Therefore, we take the transition function Ri as 


Pi=t,,-tee=2pp"[l-(P+l)'’] In 


( 12 ) 
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Taking logarithmic differentiation of eqn. (12), with respect oir, we get 

dldr(\n i?i)=|i7^+nP{ l-[p(P+l)"-'/l-(P+l)”]c?P/jp}/r ( 13 ) 

Substituting the value of dPId^ from eq. (8) in eq. (13) and taking the asymptotic 
value as 

P -> -1, after integration one obtains 

Pi='C;.;.-tee=^! |Lir'^"exp/(r), (14) 

where is a constant of integration and 

/ {r)=\[{n-\)clr-{c' l\-c)^'”'\dr,c = dcldr,c= 2\ilX +2|a. 

Substituting eqn. (14) into eqn. (7) and integrating 

x„ = Ai-Ax Ix{r), (15) 

and Tee = Xrr - ^ 1 exp /(r), (16) 

where A 2 is a constant of integration and 
/,(r)=J^r-'"-’exp/(/')Jr. 

The asymptotic value of P, obtained from solving eqn. (12) and (14) as 

p” = (^i/2) K exp/(r) (17) 

In view of eqs. (6) and (17), other stresses in the cylinder are 

Xzz={(l-c)/2-c}[ T„ + xe&]+{2d2/n-A]S3}S2/r, (18) 

and Tez = ^iiS'4, (19) 


where S 2 = {(3-2c)/ 2-c} fir, S 3 = a" r'" exp /(r), 
S 4 = exp /(r) 
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Applying bounding conditions (10) and (1 1) to eqn. (15) and (18), we obtain 

A2=AMci).Ai=^ -Idih In / 3 , ( 20 ) 


where 


a a 

h = 1^2 dr, /3 = J {[/i {a)-h{r)]rh' l{ 2 -cf -S2S, )}dr 

0 0 


The twisting couple in the cylinder is given by 


a 

M=lnA^ jr^S4dr 

0 


( 21 ) 


Eqn. (15), (16), (18) to (21) provide transitional stresses for stationary state of 
creep for a cylinder made of the material having variable modulus of rigidity, 
subjected to torsion. 


For a cylinder having non- homogeneity in accordance with eqn. (2), the 
transitional creep stresses obtained by substituting eqn. (2) into eqn. (15), (16), (1 8) to 
(21) as 


Trr = ^3[/4(«)-/4(^')]. 

T06= .^3[/4(a)-/4('-)-5'5 ], 

'tzz~-^3‘5'6[2l4(Cl)— 2l4(?' )— S’5]+iS'7(2po<5l2/ W— zfsOt^iS'g). 
tez =A3a”'^ S9, M^lnA^a"'^ h, 

where 


( ( 22 ) 


J 


A-i-2\i^d2{a^'’^ll-k+l(,ynIi,l4{r)= JiSj/r dr,l2= ‘S'9 dr, 

0 
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5'5=r''exp/,,5'6=//2/+r,5'7=(3+/r'*)/2r'^+A Ss^r'" exp/i, 

S9=r‘^ exp/i, 5',o=/^V(2r^+0^5i,=57 r'’" exp /,/=2po/^ 
ti=-k-2n,t2=-k-3n/2,fi=in-\)tldr/r(r^+t)+kt/do"kr'’-'//+t]dr. 


^-^do/r,do being a constant. 

For an incompressible material, that is, X -> eqn. (22) gives 

T:,,=A^[a'^-il+ti/2)r‘^]/ti+{Eod2/n-2A^a"r-V2)r-’‘ , 
'tez-A4a"'^r‘^,M=2‘nA4a'''^a‘^ /ti, 


I (23) 


j 


where po -> Eo/3, Eq is the constant value of Young’s modulus, 
t3=3+t2,A4=2(2-k-n)Eod2a'’ni2-k)n a”. 

For homogeneous case, that is. A: = 0, Cq is a constant and eqn. (14) reduces to 
Xrr-T:8e=Ai r‘\ (24) 


where t 4 =(«-l)Co- 2 «. 

Consequently, the transitional stresses for steady state creep become 


,(a'^-r"')/t4,X6e=^ i («"'-( 1 +?4)r'']/t4, 
x„=^i(l-Co)[2a"'-(2+t4)r"']/t4(2-Co)+{3-2co}[2po d 2 ln-Ai a" r'^]/2-c 
Xez = Ai M = 2% Ax a‘^ / h , 




°y ( 25 ) 


j 


where 


t5=(«-l)Co-w,t6=(n-l)Co-3«/2,t7=3+/6, 

.4i=(2+t5)po 
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Numerical Discussion 

In the specific case, calculations have been made for the shear stress ratio x= xe^ / 
Ml, where M\ = MU nc^, for different combinations of k and «, using eqn. (23). Fig, 1 
shows variation of the creep shear stress ratio, x with radii ratio, R = ria of the non- 
homogeneous incompressible cylinder under torsion. It can be observed that the shear 
stress is maximum at the outer radius of the cylinder for the non-homogeneity ;t<0 
with smaller values of the measure «, 0 < « < 2/3. For A:>0, the shear stress is found 
to be maximum at the centre of the cylinder and its value increases abruptly at the 
centre with increasing k and n. However, the same shear stress distribution can be 
obtained for various values of k and n, and the stress remains constant through the 



Fig. 1 Variation of shear stress in a non-homogeneous cylinder under torsion 
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cylinder for the combinations n = 1/3, A:= -1/2 and n = 2/3, A=-l . This may provide a 
choice to the designer in selecting an appropriate material for the manufacture of 
cylindrical components. 

Numerical calculations reveal that the shear stress ratio reduces to zero for the 
homogeneous material (^=0) with «=2, as the torque M required in this case, 
theoretically becomes infinitely large, and the stress appears to be negative for k>Q 
and « > 2. More realistic knowledge regarding the effects of non-homogeneity on the 
creep stresses of cylinders of various materials, subjected to torsion can be obtained 
from the numerical discussion of eqn. (22). For example, taking Poisson’s ratio Uo = 
1/3, the value of t from eqn. (22) comes out to be unity and integrals may be 
evaluated by Simpson’s rule for suitable values of k, n etc. 

Thus, the effects of non-homogeneity of the assumed form on the torsion of 
cylinder in the theory of creep are prominent. The presence of non-homogeneity A:<0, 
reduces magnitude of the shear stress significantly and may be beneficial from a 
design point of view, while for ^0, the stress increases promptly at the centre of the 
incompressible cylinder and may not be useful with higher values of k and n. 
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Abstract 

The dynamic problem in micropolar viscoelastic medium has been investigated by employing 
eigenvalue approach after applying Laplace and Hankel transformations An example of 
infinite space with concentrated force at the origin has been presented to illustrate the 
application of the approach The integral transforms have been inverted by using a numerical 
technique to obtain the radial displacement, normal displacement, normal force stress, shear 
stress, couple stress and microrotation in the physical domain The results for these quantities 
are given and illustrated graphically. 

(Keywords micropolar viscoelastic medium/eigen value/Hankel and Laplace transformation) 


Introduction 

Modem engineering structures are often made up of materials possessing an 
internal structure. Polycrystalline materials, materials with fibrous or coarse grain 
structures come in this category. Classical elasticity is inadequate to represent the 
behavour of such materials. The analysis of such materials require incorporating the 
theory of oriented media. For this reason, micropolar theories were developed by 
Eringen'’^ for elastic solid, fluids and further for non-local polar fields and are now 
universally accepted. A micropolar continuum is a collection of interconnected 
particles in form of small rigid bodies undergoing both translational and rotational 
motions. Different authors'*"'^ discussed different type of problems in micropolar 
elastic medium. 

Eringen'^ extended the theory of micropolar elasticity to obtain linear constitutive 
theory for micropolar material possessing internal friction. A problem on micropolar 
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viscoelastic waves has been discussed by McCarthy and Enngen'^. They discussed the 
propagation conditions and growth equations governing the propagation of waves in 
micropolar viscoelastic medium. Cicco and Nappa'^ discussed a problem of Saint 
Venanfs principles for micropolar viscoelastic bodies. Kumar et al}'' studied Lamb's 
plane problem in a micropolar viscoelastic half-space with stretch. Recently, Kumar'* 
discussed wave propagation in micropolar viscoelastic generalized thermoelastic 
solid. Most of the problems studied so far, in micropolar viscoelastic, involve the use 
of potential functions However, the use of the eigen value approach has the advantage 
of finding the solution of eqns in the coupled form directly, in the matrix notations, 
whereas the potential function approach requires decoupling of eqns. Yet, much work 
has not been done in micropolar viscoelastic solid using the eigen value approach. 
Mahalanabis and Manna'® applied eigen-value approach of linear micropolar elasticity 
by arranging basic eqns of linear micropolar viscoelastic m the form of matrix 
differential eqn. Recently, Mahalanabis and Manna*® discussed the problem of linear 
micropolar thermoelasticiy by using eigen value approach 

In this paper, we consider a two dimensional axisymmetric problem in a 
homogeneous isotropic micropolar viscoelastic medium. The solutions were obtained 
by using eigen value approach after employing integral transformations technique. 
The integral transforms were inverted using a numerical approach. 

Basic Equations 

Following Eringen'"* the constitutive relations and the field eqns in micropolar 
viscoelastic solid without body forces and body couples can be written as 


t,j - XiUkx^,j + + K, {uj^,~B,jk<9k), (1) 

m,j = aiijit.i h,j + pi y, (2) 

(X, +2|a, +A’,)AA.m-(p, -i-^r,)Ax Axm-i-.K,Ax$ = p — — (3) 


(a, + P, + y, ) AV$ A X M - 2^,$ = p; ^ 

or 


( 4 ) 
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where, 

X, = L + = ^ + Kx=K + Kv— (5) 

dt dt dt 

a!=a + ar— , P!=P+PkA, y, = y + A (6) 

dt dt dt 

where, L, p., ^ , a, p, y, Xv, IIk. K-vO-v, Pk, Yr are material constants, p is the density, j is 
the micro inertia, u is the displacement vector, | is the microrotation vector, t,j is the 
force stress tensor, m,j is the couple stress tensor. The necessary and sufficient 
conditions for the internal energy to be non-negative, as given by Eringen''* are 


0 < 3 A, + 2p. + K, 

0<|a, 

0^^ 


0 < 3a +2y 

-y ^ p < y, 

0 ^y 

(7) 


and 


0 < 3 Xy+ 2|Xv 

0<p, 

0<K 


0<3av+2jv 

-jv^^v^yv. 

0 <yv 

(7) 


Formulation and Solution 

We consider a homogeneous, isotropic micropolar viscoelastic medium of infinite 
extent with cylindrical polar co-ordinate system (r,9,z). To analyze displacement and 
stresses at the interior of the medium due to concentrated load, the continuum is 
divided into two half-space defined by 
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»Half-space I: r <oo 


-00 <z < 0, 


»Half-space !!:/•<<» 


0 ^ z < 00 . 


As we are discussing an axi-symmetric problem so we take 


U = {Ur, 0, Mj), 


(0, (l)e, 0) 


( 9 ) 


We define dimensionless quantities as 




z 

v 


, ph(o'\ 
u. =- - 


1^1 


, phca'^u, 
u=- ‘ 




<Pe = 

Pi 


pAV ’ 


t’ =—t , 

a ^ 


t' = —t_ , 

a 




K, 




CO 


Pj 


( 10 ) 


where A is a constant having dimension of length. 


Due to axi symmetry about z-axis the quantities are independent of 9. With these 
considerations and using (9) and (10) in eqns. (3) and (4) we obtain : (on suppressing 
the primes). 


1 a 

dr^ + r dr 


u, + 


(!-<,=) 


d 2 ^ 

+ J. 


drdz 


dz 


2 *^4 




1 d\. 


dz (^i+j'j) St 


( 11 ) 
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^ \drdz r dz \ \_dr^ r dr \ V (5, +^2)5^^ ’ 

( 12 ) 


g^(j)a 1 5(|)e ^^<|>e (|)e 2h^s^ ( du^ ^z'] _ 


_JLH_J LE.J L£. — __ , 3 --r -■'■'2 I Te /i7\ 

‘I a-' r dr &= rM J j Uz ITJ'!?" ' ' 


where 


p/t"(X2 +11/) CO* 




ph^ K, co*^ 


yiPh^ CD* 


*4=^-T 


2 *^2 
I = ^ 

S,+Jj 


5, +5j 


Applying Laplace transform w.r.t. time variable Y defined by 
K ('•. P)> “z if‘> P% • 2 . P) = 

00 

J{m, (r, 2 , 0, Wz(r, z, 0, <t>e (?•, 2 , dt (14) 

0 

and then Hankel transform w.r.t. space variable V defined by 

'^^{^,z,p)\='^\uXr,z,p){z,r,p),u,{r,z,p), $9 {r,z, p)]rJo{%r)dr, 


00 

= J“r (r,z,p)rJi (^r)dr, 


( 16 ) 
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on eqns ( 1 1 ) - (13), we obtain 
1 


r 

\ yl , P 

- .1Li^ 

s + 

Mr + 2 

L 

a 


(17) 


u, =\a c, + 

‘ [ 51+^2 




(18) 


and 


Kjh^ Kih\ I 

= i Uz + 

Y/ Y/ 


2 

s2 4. 1 + £_ 

Y/ ^4 


<(>0 » 


The system of eqns. (17) - (19) can be written as 

— Wi^, z,p) = A ii„p) W {^,z,p), 
dz 


(19) 


(20) 


where, 



'ul 


'0 /] 


Mr 

W = 

[v 

A = 

5 



’ u = 

“z 






_4’9. 


0 = 


A = 


-(l-a^k 


a' 


Yi 


0 

0 



0 



'1 

0 

o' 

■q 

0 


0 

0 

0 

1 = 

5 


1 

0 

_o 

0 

0_ 


L» 

0 

1 
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A = 


a 


2^ P 


2 \ 




(^1 +^2). 


0 

0 


To solve eqn. (20), we take 

W(l^,z,p)=Xa,p) 

so that 


2e , P 


a% + 


(5, + ^2 ) 


ml 

Yi 


.2 , 2Ky ^ p^ 

Yi -^4 


( 21 ) 


( 22 ) 


A(%,p)W(^,z,p) = qWa,z,p) 


(23) 


which leads to eigen value problem. The characteristic eqn. corresponding to the 
matrix^ is given by 

det [A-ql] = 0 (24) 

which on expansion provides us 

q^ - Xiq^ + Xiq^ -Xi = Q 


where 




1 + 


o' J 




2K,h^ 


(i’l +52)+[ y, 


J 


^2 


2K^^p^ 
Yl ^4 . 



1 + -1 
V a 



(26) 
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Y.a L 


-7-^ — \ +— — \ — \ 

(s.+ijJ (‘yi+52)jL (^i+'^a) 


1 .2 , r 


a' ^ {Si+S^)j L (^,+^2)J L Vi 


2 . 2K,e . p^ 


sXh%\ 


(s,+5j 


The roots of eqn. (25) are ± 9 ,, /= 1, 2, 3. 

The eigen value of the matrix A are the roots of eqn. (25). We assume that real 
part of g, are positive. The vector X (^,p) corresponding to the eigen values q, can be 
determined by solving the homogeneous eqn. 

[A-qI]X{%,p) = 0 (29) 

The set of eigen vectors AT, i^,p), {i= 2, 3, 4, 5, 6 ) may be obtained as 


(k,p) 


tp) 

tp) 


where, 


^.ia,p)= b, , 


X.2 i^,p)- 


a.q. 

^q, 


q=q„ /= 1, 2, 3(31) 
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1 

fN .... 

I 



A), = 

b, 

. J 

, 2^2 {^,p) = 

1 

J 

1 

, y=z+3, q=-?„ /= 1 , 2,3 

(32) 






Yi 


2^* 

4 


"4 J 


I Y, 


(33) 


b. = 


2 \ / 


( 

(^i + S 2) 

r „2 ^ 

(•^1+^2). 


,2 

4 + — - — + — 
Yi -54 


>1 1 

N f 

+ a 

X, - 

J 1 

l V 


.2 , , P' 


Yi 


^2 , K,s\h\] 

Hi •“ 

Yi 


(34) 


A = 




Y, 


■■ 

2 M 



^ (5j+52)J_ 


(35) 


The solution of eqn. (20) is given by 

W{%,z,p) = i [B,X, (^,p) exp (g,z) + (^,p) exp (-^„z) 


where, B, (z = 1, 2, 3, 4, 5, 6) are arbitrary constants. 


The eqn. (36) represents the solution of the general problem in the axisymmetric 
case of homogeneous isotropic, micropolar viscoelasticity by employing the eigen 
value approach and therefore can be applied to a broad class of problem in the 
domains of Laplace and Hankie transforms. 
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Application 

We consider an infinite micropolar viscoelastic space in which a concentrated 

force of magnitude F = acting in the direction of the z-axis at the origin 

° litr 

of the cylindrical co-ordinate system The problem is axismmetric w r.t. the z-axis. 
The boundary conditions at the plane z = 0 are given by 

M, (r,0", 0-«r (r,0',O = 0, Wj (r,0\ 0-Wr = 0, (37) 

(1)9 (r,0", 0- <|)6 {^0',r) = 0, (38) 

t)- m-e {r,0',t) = 0, (39) 

tzr{r,0\ t)- tzr (r,0-,t) = 0, (40) 


Uz(?'>0^ 0- tzz ir,Q',t) 


Fo6(r)8(t) 


Applying the Laplace and Hankel transforms to eqn. (37)-(41), 


we get 


-M,(^,0» = 0, M,(^,0»- u, i^,0-,p)=0, (42) 


?e(^.0\;.) - (j)e (^,0',/?) = 0, 


,P) - w^(^.0',;7) = 0, 
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The transformed displacements, microrotation, couple stress and stresses are 
given for z > 0 by 

("^1^) (““92^) ^ 3 ? 3'86 CXp ( 47 ) 

{^^,p) = biBi exp (-^12) +62^5 exp (-922) + ^>3^6 exp {-932) ( 48 ) 

= ^[-84 exp (-^iz)+ 55 exp (-92z)+56 exp {-qzz)], ( 49 ) 

(^,z,;7) = exp (rqiz)+q2B2 exp (-q2z)+qiB6 exp (-qiZ)], ( 50 ) 

L = («' "'s -48 i 59 + ^!o) 84 exp H,r) 

+ {a2 qls^ -^b2S9 + ^10) -85 exp {-q2z) 

+ (<33 ^3^8 -^8359 + ^■sio) Bi exp i-q^z) (^ ^ ) 

^ (^^,J3) = -[91 ( 8 i 56 + ^ai^s) Bi exp (-qiz) + q2 {b2S6 + ^^255) -85 exp (-^2^) 

+ 93 (83S6 + ^<3355) -86 exp i-qsz), (^^) 

and for z < 0 by 

(^^,/j) = 8 i 5 iexp (?iz) +82^2 exp (^22) + 83^3 exp (93^) (^ 4 ) 
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= -^[B\ exp {q\2)+B2 exp (g2z)+-63 exp (qiz)], 

= -siiXq\B\ exp {q\z)+q2Bi exp {q2z)+qzB2 exp (932)], 


?; {%^,p) = (a, q^,s, + ^^o) 5, exp (q^z) + (fl2 qls, -^b 2 S 9 + ^s,o) B2 exp {q^z) + 

(03 qls^-^b^sg + ^ 5 io) Bi exp (q^z) 


T„ {^^,p) = q\ {b\S(, + ^aisi) B\ exp {qxz) + qi (^>256 + '^aiSs) Bi exp {q2z) + ^3 (& 3^6 + 
^<3355) 53 exp {qiz), 


where, 


= . ^ (^i + 2 p, +i:i)pi S 7 = __li 


S% 


^ (p,+i:,)p, 

pA:,/z'(o‘' 






2 , *2 ’ 


5 io 


= 




p/z 0) 


Using conditions ( 42 ) -( 46 ) in eqns. ( 47 ) - ( 48 ), we obtain 
a\q\ {B\ + ^ 4 ) + Qiqi {Bi + 265 ) + <^3(^3 (-S3 + B^ = 0 , 


{Bx - S4) + (52 - 55) + (53 - 56 ) = 0 , 


b\ ( 5 i “ 54) + Z?2 {Bi - Bs) + bi (53 - 55) - 0 


qx (5i + 54 ) + qj (52 + Bs) + qi {Bs + 56) - 0 

[a, g-f Sj -^bxS9 + ^10] (j?i - 54) + [fl2 qls^ -^b 2 S 9 + ^Sio] {B2 - Bs) 


( 59 ) 

( 60 ) 
( 61 ) 
( 62 ) 
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+ (<33 + ^lo) (53 - Be) = 0 


(63) 


q\ {b\S 6 + ( 5 i + 54) + q2 {biSe + ^0255) {Bj + Bs) 


+ qs {bsSe + (B^ + Be) 

2 % 


(64) 


Solving system of eqns. (59) - (60), we obtain 


B\ - B 4 


^0 (^3 ^2 ) 

47tq',A, 


(65) 


Be — B$ 


Bi - Be 


4tC^2^1 

471^3 A, 


( 66 ) 


(67) 


where 


A] = S 6 [( 02 ^ 3 - a^bi) + (<33^1- ^1^3) + (« 1 ^ 2 - Oib]). 


Thus the components u^, w, , 4»e > ’^nd 7^ have been determined in the 

transform domain and these enable us to find the displacements, microrotation, couple 
stress and forces stress. 


Inversion of Transforms 

The transformed displacements and stresses are functions of z, the parameters of 
Laplace and Hankel transforms p and ^ respectively, and hence are of the iormfi^^p) 
in the physical domain, first we invert the Hankel transform using 
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0 

Thus, expression (68) gives us the Laplace f (f^,p) of the function We 

invert the Laplace transform by using the inversion technique as used in Kumar and 
Choudhary^' 

Numerical Results and Discussion 

Following Gauthier^^ we take the following values of relevant parameters for the 
case of aluminum epoxy composite as 

p = 2.19 g/cm^ A, = 7.59 x 10’“ dyne/cm^ 

p = 1.89 X 10‘° dyne/cm^ K= 0.0149 x 10’“ dyne/cm^ 

Y = 0.0268 X 10'“ dyne, j = .00196 cm^. 

For a particular model of micropolar viscoelastic solid the relevant parameters are 
expressed as 


X = X (l + *2r') / = 1, 2, 3, 4 for X = A, p, K, Y respectively 

where, 

01 = 0.05, 02 = 0.01 03 = 0.015 04 = 0.1. 

The comparison of values of dimensionless normal displacement Uz [=(27r«z/Fo] 
and couple stress Mg, [=(2tc AWiv] for micropolar elastic solid (ME) and micropolar 
viscoelastic solid (MVE) have been studied. The computations were carried out for 
two values of dimensionless time 0.05 and 0.10 and for A = 1 at z = 1.0 in the range 0 
^ r ^ 10. The solid lines [ — ] in graphs represent the variations for time 0.05 and the 
dashed [ — ] represent the variations for time 0.10. For both the times, solid or dashed 
lines with center symbols represents the variations for ME, whereas without center 
symbols are for MVE. 
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Fig. 1 shows the variations of normal displacement Uz with r. In the range 0 < 

1 and 5.5 <r< 10 the values of for time 0.50 are less in comparison with those for 
time 0.10 whereas in the range 1 < r < 5.5 the time 0.05 are more than for time 0 10 
for both MVE and ME. For fixed time 0.05 the values of for MVE are more than 
those for ME in ranges 0 < r < 1 and 5.5 < r < 10 and less in the range 1 < r ^ 5.5. 
The similar trends for time 0.10 have been shown. 

Fig. 2 shows the variations of couple stress with r. The values of A4e have 
been taken by multiplying their original values with 10^. In the ranges (0 < r < 2 ) and 
7 < /• < 9 the values of for time 0.05 are less than the values for time 0.10 and in 
the ranges 2 < r < 7 and 9 < r < 10 the values for time 0.05 are more than those for 
time 0.10 for both MVE and ME. For fixed time 0.05 the values of M^e for MVE arre 
more than those for ME in the ranges 0 < r < 2 and (7 < r < 1 0) and less for ME in the 
range 2 < r < 7. The similar trends for time 0.10 have been shown. 
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Abstract 

The Kelvm-Helmholtz instability of two superposed viscous conducting fluids in a uniform 
vertical magnetic field is discussed m the presence of effects of surface tension and 
permeability of porous medium. The dispersion relation has been derived through normal mode 
technique and solved numerically The streaming velocity is found to have a destabilizing 
influence while viscosity, surface tension and permeability of the medium have all stabilizing 
influence on the growth rate of unstable mode of disturbance 

(Keywords Kelvin-Helmholtz mstability/surface tension/ viscosity/porous medium) 

Introduction 

The problem of the Kelvm-Helmholtz discontinuity between two superposed 
fluids is of prime importance in various astrophysical, geophysical and laboratory 
situations. The Kelvin-Helmholtz discontinuity arises when air is blown over mercury 
or when highly ionized hot plasma is surrounded by a slightly cold gas or when a 
meteor enters the earth's atmosphere. Chandrasekhar^ has given a detailed account of 
the various investigations of this problem as investigated by different researchers for 
incompressible fluids. For the case of a low - P magnetized plasma m the presence of 
effects of compressibility, D' Angelo^ found that a velocity difference between 
adjacent field-aligned flows, on the order^of the ion-acoustic speed Cg, is sufficient for 
instability. The results of his calculations were in. agreement with experimental 
findings in a Q-machine (D' Angelo and von Geoler^). Gerwin"^ has examined the 
stability problem of the non-conducting streaming gas flowing over an incompressible 
conducting fluid. Kalra, Lister and Hosking^ have investigated Kelvin-Helmholtz 
instability of superposed fluids in the presence of the effects of magnetic resistivity 
while Bhatia and Steiner^ have analysed this problem for partially ionized plasmas. 
Shivamoggi^ has studied the effect of finite magnetic resistivity on the Kelvin- 
Helmholtz instability of superposed plasmas. Sengar^ has analysed the stability of 
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superposed gravitating streams in a uniform vertical magnetic field in the presence of 
effects of magnetic resistivity. Mehta and Bhatia investigated ICelvin-Helmholtz 
instability of two viscous, superposed, rotating and conducting fluids. Singh and 
Khare'° have investigated the stability of two semi-infinite homogeneous gravitating 
streams of infinite conductivity under uniform horizontal magnetic field and uniform 
rotation. Eldabe and Hassan" have examined the effects of surface tension and finite 
magnetic resistivity on the Kelvin-Helmholtz instability of superposed fluids in a 
horizontal magnetic field. Bhatia and Hazarika'^ discussed this stability problem for 
superposed viscous gravitating fluids The importance of the Kelvin-Helmholtz 
problem has been demonstrated recently by Benjamin and Bridges'^"’'' who have given 
an excellent reappraisal of the classic Kelvin-Helmholtz problem in hydrodynamics. 
They have shown that the problem admits of a canonical Hamiltonian formulation and 
obtained several new results. More recently Sharma and Kumar'^ have studied the 
Rayleigh-Taylor instability of two superposed conducting Walter's B' elesticoviscous 
fluids in hydromagnetics while Allah'® has investigated the effects of magnetic field 
and heat and mass transfer on the Kelvin-Helmholtz instability of superposed fluids. 

The flow through porous medium has been of considerable importance in recent 
years particularly among geophysical fluid dynamicists and petroleum engineers. The 
physical properties of the comets, meteorites and interplanetaiy dust strongly suggest 
the significance of the effects of porosity in astrophysical context (Me Donnel'’). 
Several authors'®'^^ have investigated the effects of permeability of the porous 
medium on different stability problems. 

The aim of this paper is to study the Kelvin-Helmholtz discontinuity between two 
viscous conducting fluids in a uniform vertical magnetic rield through a porpus 
medium in the presence of the effects of surface tension. 

Perturbation Equations 

Consider the motion of an incompressible, infinitely conducting viscous fluid of 
variable viscosity \xo(z) through porous medium having a streaming velocity U along 
x-direction. The magnetic field is taken to be uniform and vertical. If we denote by 
u (u, v,_w), h {hx, hy, hz ), 5p and bp, the perturbations in Darcian velocity, magnetic 
field H , density p and pressure p respectively, then the linearized perturbation 
equations appropriate to the problem are 

P ^ + £(Lf.v)M =-W5p + g5p + (vxh)xH+ 

e ^ .r 4, r V / g 
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6— +(c7.v)A = (^.vk 

dt ^ ^ ^ ^ 

(2) 

e -^Sp + (t/.v)5p + = (w.V).p = 0 

(3) 


(4) 

V. w = 0 

(5) 

V.^ = 0 

(6) 

As stated above, we take streaming velocity U = (17,0,0) and magnetic field H = 
(0,0, Ho). Here, T, p and g = (0,0, -g) denote surface tension, coefficient of viscosity 
and acceleration due to gravity respectively. In equation (1), 5(z - Zs) denotes Dirac's 5 
function, X denotes permeability of porous medium and e is the medium porosity. 

Analyzing the disturbances in terms of normal modes, we assume that the 
perturbed quantities have dependence on space coordinates x, z and time t of the form 

i^(z) exp (ik,cX + nO 

(7) 

where F(z) is some function of z, is the wave number along x- 
be complex) is the frequency of the harmonic disturbance. 

-direction and n (may 

Using equation (7) in equations (1) to (6), we get 


+ ikU) u = -ik„hp + H,Dh^ - HJkX + -(^' 
e e 

-k^)u--u (8) 
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— (« + ikp ' ) V = H^Dh^ + — {d^ - kl 

e ^ 

)v--^v (9) 

{n + ik^U)w = -D?)p-g&p + ^{D^ -1 

w--^^6(z-zjw (10) 
A, ne 

e (n + ikx U) hx = HqDu 

(11) 

€ (« + ikx U) hy = HqDv 

(12) 

G (« + ikx U) = HqDw 

(13) 

(n + ikx U) 5p= -(Dp) w 

(14) 

ikx hx + Dhi - 0 

(15) 

ikx M + Dw = 0 

(16) 


where 



Eliminating some of the variables from the above equations, we obtain an 
equation in w as 

—[klpw - D{pI>A-^w + ^(d^- kl 
€ n' e n'e 

+^iD^-kl)\-HD^-kl)w-v^b{z-z )w = 0 (17) 

€ A, n'e 


where we have written «' = (« + ikJJ). 
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Two Superposed Streaming Fluids 

We now consider the case when two superposed fluids of uniform viscosities |lii 
and |i 2 ? uniform densities pi and p 2 , magnetic fields H\ and Hi and with streaming 
velocities U\ and Ui are separated by a horizontal boundary at the interface 2 = 0 
Therefore, in both the regions 2 < 0 and 2 > 0 of constant densities, equation (17) 
becomes 


(D^- kl) -M^) w = 0 


(18) 


where 


M‘ 


'2 _ 


kln'v 


«'v + F? 


1 + 


6 n 
Xkl vk^ 




(19) 


IX H 

In equation (19), v = — is the coefficient of kinematice viscosity while Vj = — ^ is 

P P 

the Alfven velocity 

Now, seeking the solutions of equation (18), which remain bounded in the two 
regions, we obtain 


w, = Ax n[ e'^+Bi n\ ; z < 0 (20) 

W2 = A2n[ +B 2 n[ ;z>0 (21) 

where Ax, A 2 , Bx and B 2 are constants of integration and Mx and M 2 are square roots of 
equation (19) and 

n[ = n + ikxUi i n[ = n + ikx U2 ( 22 ) 


and 
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M^ = 


kln[v 


Mi = 



1 + — ^ 

n[v,+V,r 

X^k. 

kln\vi 

1 + -^ 


-+- 




«2V2 + 


X^ki Vjk 


2’^x 


(23) 


(24) 


are the values of Mi and Mz for the two regions. Since «’ is complex, Mi is therefore 
complex and consequently M\ and Mi are complex. In writing the solutions (20) and 
(21) it is assumed that M and Mi are so defined that their real parts are positive. 

The above solutions must satisfy certain boundary conditions. These conditions 
require that at an interface z = 0, 

w, Dw and p. (Z)^ + kl)w (25a, b c) 


must be continuous. 

If we integrate equation (17) across the interface z = 0, we obtain another 
condition 


n'lP^Dwi - p2 [d^ - kl )Dw^ ) ^2 

n-, A 


J2=0 


nlp^Dwi -k^x)Dw^ -^(d^ 

A 


J2=0 




\K «i'y 


-klT\ wo-2kliiii-p,)iDw)o 


(26) 


where wo and {Dw)o are the unique values of these quantities at z = 0. If we apply the 
conditions (25) and (26) to the solutions (20) and (21), we obtain 


Ai + B\ = Ai + Bi 


(27) 
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kx Ai + MiBt = -kxA 2 - M2B2 

^2jr2 


2 nj 


■kl (ia .2 -1^1) (Ml + MiBx - M2 - M2B2) 


Eliminating the constants Ai,B\,A 2 and B 2 from equations (27) to (30) and writing 


SK ( ^2 «. 1 Kt 


2 U 2 ^(P.+Pz) 


C = klia2V2 -aivi), ai.2 = 


P.+P 2 


we obtain the following characteristic equation 


(Mi-kx) 2kl{ay, -a,v 


+ - kl ) {a,n; +a,n[- 2 i?}] 



504 


P.K. BHATIA and ARUNIMA SHARMA 


a,v, (M,^ + — (Mj -A:^) a^n' 


+j{M, -K)- a,n[ 


-2k, a,v, {m^ +-^(^2 -^x)-« 2«2 


+ V , {Ml -kl)\ ^ M, + ^ (M, - A:J - a,^' 


+{Mi-k,) [a,v, [mi - kl ) {a + a,«; - 2R} 


-2A*(a,v,-a2V2 =0 (33) 

I «. K JJ 

wherein the subscript A on Alfven velocity Va has been dropped for convenience. The 
dispersion relation (33) is quite complex. We therefore, consider the case of two 
highly viscous and highly permeable fluids for then we can Avrite 


aj«j +ajV, 


M, = A + 


< 


2k a,v, + — L- 
‘ ' 2n[ 


a2n2 +a2V2 


M2 = k + 


62 k^V 2 
X2 n'2 


2k a,v, +-^ 
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Now we consider the case of two streams, of equal Alfven velocities and same 
kinematic viscosities and same permeabilities, flowing past each in opposite directions 
with equal velocities i.e. 

Fi = F 2 = V, V, = V2 = V, Ui = U= U 2 = -U 

ei = €2 = e, A-i = ^2 = ^ (36) 

Now substituting the values of Mi and M 2 from equations (34) and (35), the 
dispersion relation (33) becomes, on writing kx = k, the ninth degree polynomial 

(37) 

/=0 


where the coefficients A, (i = 0, to 9) are given in the Appendix. 

Discussion 

The dispersion relation given by equation (37) is quite complicated particularly as 
the coefficients A's are complex. We have therefore solved the dispersion relation 
numerically for different values of the parameters characterizing viscosity, streaming 
motion, surface tension and permeability of the medium. The dispersion relation 
was of course first non-dimensionalized by measuring n, v, B, V, Q and U in terms of 
'ig- For numerical evaluation of the roots of dispersion relation we have taken 
tti = 0.25, a 2 = 0.75 (the unstable configuration), Alfven velocity F = 1 in all the 
calculations. 

These calculations are presented in Figs. 1 to 4 where the growth rate (positive 
real part of n) is plotted against the wave number k for several values of v, U, B and 
Q. In Fig. 1, the growth rate is given against the wave number for v (viscosity) ~ 0.1, 
0.2, 0.3 taking fixed values of B (surface tension) = 0.1, U (streaming velocity) = 1,Q 
(permeability) = 0.05. From the Fig. 1 we find that the growth rate decreases on 
increasing the viscosity for same values of the wave number. We thus find that 
viscosity has stabilizing influence on the Kelvin-Helmholtz instability of the 
superposed fluids. In Fig. 2 we have given the variation of the growth rate against the 
wave number for the values of streaming velocity U = 2.1, 2.2, 2.3, 2.4. In these 
curves also the values of the other parameters have been kept fixed. From Fig. 2 we 
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find that the streaming velocity has a destabilizing influence as the growth rate 
increases with increasing streaming velocity for same values of the wave number. 
Fig. 3 depicts the variation of the growth rate with surface tension. The curves in the 
Fig. 3 clearly show that the growth rate decreases with increasing values of B, the 
parameter characterizing effect of the surface tension. As the growth rate decreases on 
increasing the values of B for same wave number k. The effect of surface tension is 
thus stabilizing on the growth rate of the unstable configuration. Fig. 4 gives the plot 
of growth rate against the wave number for the values of permeability parameter Q = 
0.2, 0.1, 0.05. From Fig. 4, we find that the grovrth rate decreases on decreasing 
parameter Q. The effect of permeability is consequently stabilizing and/or the effect of 
porosity is destabilizing. These observations are in agreement with the earlier 
observations of several researchers oh the different stability problems. 



Fig. 1-Variation of growth rate n against wave number K for v= 0.1,0.2,0.3 
taking 5=0.1, iA=1.0, g=0.2 
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Fig 4~Vanation of growth rate n against wave number KfotQ- 0.2, 0 1, 0 05 

taking U=l.0,B-0 1 
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2 ia^a 2 - 9 UX(aj,-a^) + 

+ 2a^aj, 2-9+ 27^11 - 2a^a^)] 


+ 2 aj^lc^) 


= 2 iUk 9 a [ 2 “i ®2 ^ 

+ 4 o|a| + 0 if 2 {^ f 2 Q(Q+ 2 k‘^C 3 a ^-a^) 

+ 4 a,a,^VBQ + ^ V [4Q( l-cc ,a,) - k^4a .a^iS- 12a 

+ ag-a^j-e)] + 

Ag = iUkaj^agl" ^ 9^k2u2a^a2U2-a^)+2a^a2^^QL2k^Ua^-a2) 
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+ ^^Q[V^i3a^a2Q-4k2u2-a^)ll-2a^a2)) 

- 2aj^a2kCa2“a^)+2aj^a2^ Qk^^kB-ii^ )]+ [kCsg-a^) 
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- Y^-[lc^BV^(l-2aja2)+-i)^a^a|kQ(3Q-4k^)+ 9V k [k 

. (2-9a j^ci2+4a^a2)+ a j^a2'4]+ 

- a tag-a ^1 + ^ k V { k^ [{<x^-a ^) (7a ^ag-^ ) 

+ 4aia2(l-4a^a2)-2]+ (aa-a [4k^3a ^ag- l)-4a ^ag+Q]} 

+ 4a ^ag t) Vu^+a ^agkV { -25 9 [Q(2-5a ^ttg) 

- 2k^( 1 - 20 ^ 02 )]+ ll[v'^+-i)^V^[Q(a^a2+2)+4k (a^ttg 

. (a2-a^-12a^a2+7)-l)> 2^ k2B(2a ^ag -9 Vv^( 1-a ^ 

+ 9^Q(aiagQ+2a^agk2(3a^-a2))]} +2kV-[ a^^ag { [-PBk^ 

. ( 2 Q- 3 k 2 )- 2 Y\^-k(ag-a^)] + ^ \q f 4 Y^k [k^ (Og-a ^-Ea^ag) . 

+ Q( (ttg-a^) (a j^ag-l)-a^ag)]- 2a j^agQtttg-a^)! 

+ 9V{k(a2-a^)(3Qa^ag +2(k^-2Q))+ 2^k^BQ(Q 
- k^(2a^a2+l))+ V^[Q^+4k^a^( l-2a2)-k^Q(2(a^a2+a2-a 
+ 5)]]] - Y^^k^Y^(2a^ag-l)-T?V^k^BU-5aj^a2) 

+ 9^ •[aia|kQ(a2-ai)+v\^[Q(5a^a2-2a^a2-2) 

+ 2k^[(ag-a^)(l-5a^a2)+ a^a2(2a j^a2-7) + l]]]|] 

+ 2a3^cigkV( 9^Q(l-2aj^cc2)+agV^)+v‘^{9\[k^QUaj^a2 

. (3ag+a2^-a j^ag)-2)- a j^a 2 (a 2 -aj^) (Q^+2k‘^)] + 
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+ 7^ {v^k^+lc^t2aj^a2Ua2+a^)-l)-Vi!:^B[;Q-2k2-*-2aj^a2 
. (3k^-2Q)]+ 2 '?\^o^[2k^Q(.2aj^(2a2-(x^)-l) 

+ 4a^a2k^ + OgQl l-Q)-ik^( l+4a|)]]} 

Ag = iUk{ 4a^a| S)^k®U®ta2-a^)+a^a2k'^U^{^2v2u2-a^) 

. [2k^-Q(2-a^a2)] + l5['J^V^k^[2(a2-a^)t3a^a2-2) 

+ 2(2-90^02^3'^ 2 -P^7^Q[a^a2-l ] +2a ^02 S>^Q [2^^^ (2a ^-ocg) 

- Q(a2-o^)]- 20^02 -i}^k - 7*^ ( 02-0 2^)-2k^B t) (2a^a2 )1 

+ 2k^U^ { 2a^a^ ^2k^7^ (a2-a^)-f7^ f ^ [( 02 - 0 (3k^^Q-k^) 

- + 2k^(4aj^a2-l)(2(Q+k^(k^-l))+k^Q)]- k 

-9 k^BU2-a;^)| + 9^v\[2k(2a^a2-l)-3a^a2Q 
-■9 k^QB(Q-k^)(a2-a^)j} + 7^ [ -a^a| ^k (3Q-4k2 ) 

+ 7^(a2-a^) [7^k^(2a^a2-l) + ^k^B(2-a^a2)+-9^k2 
. [Q(5a^a2-2)+ 4k^ ( l+2a ^a 2 ( l-2a ^ 02 ) )+2 ( Q+Ba^Og)]!]^ 

+ 3[2aj^a2v\^(^^(a^a2-a2+Oj^)‘‘‘“i^^)]- v'^ ^ 7^“^ (og-a 

- [(“2““ 1^ ^ l-a^a 2 )-a ^ k 2 B(a 2 -a (Q- 2 k 2 )] 

+ ^ k [a^OglQ-k^) ^y- 2 k^)-k^' 4 (, 1-20^02) (2 (Og-a ^)+a ^Og )] 

+ 27\^[(a2-a^) [(2k^+Q}(k^-a^(i2(2k^+Q))-2a^a2k'^] + 
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- 2kVa^a29{ ll[2a W U-5a )]- 2a^a2[2-j2 

. U^4<x/)+2ic3B->)+372j] +kV{5k2 9Vn2-7a^a2^2a^^a|) 

_ a^ttg-v)^ kU2-a^)+V^U--3k^)+2^ k BUU-a^ag) 

- k^) +i)^{Q^U-5aj^tig)+ 2Qk^C2U2-“i)^“l“2“^^~^“l“2^^^ 

- 8k\j^]} - 2a jag kQ(.a2-“i)l 

+ 2ajag9Qk^kB-2>)) + -?v2[:kU2-«i)^«l«2“2)^ 

+ (2k2aj(2-3aj)+Q)+29k3BQ(l.a^a2)]}] ^k^U^ 

. IBa^ag-J v\^(a 2 -a^)+ 2a j^cig-? V^ ^ V^k [tQ-k )(.(a 2 -aj) 

+ 9k^B(Q-2k^))-V^k(2Q-3k^)] + -? ^ [ kQ(Q-k^) (ttg-a j) (2-3a ^ag) 
+ 2V^^ [Q^ ( 2a j^ag-Og+a j) +2k^ (4a ^a^-a^+a j^) +3k^Q (ag-a j) 

+ 2(4a^a2“^^^‘^"^^^^5il k^V® { (Og-a j) [(y-2k2 ) 

. (l-ajag)+ aj^agCQ-k^)]^ k^B ^2^ V^k [2a ^ag(2k2+Q) 

+ (a2-aj^)(k2+Q)]+a^a2-9 [(Q-k^) (ag-a ^)- (Q-2k2)]} 

Aj^ = ikU ^8ajag-J>k®U®(a2-aj^) [Baj^ag •9^Q+V^(2-a ja2)]+ k'^U^ 

. ^ -4k^S) v'^(2-3aj^ag)+2aja2'P { V^(a2-aj^)(3k^-Q) - 
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- fQ^Cog-a^) (,a j^ci 2 - 2 )+ 2 k^Q(, tttg-a ) 

. (l-2a 2 ) -4a ^ ) - lea^agk^] -5a ^a k^B ta g-a ^ ) 

. {Q-k^)}}j - 4k^V^ ^2a^a^^7\^(l-a^a^)^V^lk^3v'^ 

. ( l-2a ^a 2 ) + |kQ ( 3 Q-4k^ ) + S) [k^ ( 2-9 a ^^a g+^a ^a^ ) 

+ a^a2Q]+ [2QC l-3a ^agj-a ^ag Ittg-a ^)] 

+ -D k^V^ k^ (^Cag-a^) (7a j^ag-2)+4a l-4a j^ag)-2] 

+ (ag-a^) [4k^(3a^a2-l)-4a^a2+Q]] - 2a^ag 

. [k^(2(a2-a^)-12a^a2+ll)+ Q(a2-a^)] Jj +2 [a ^^ag-^) V®k^ 
. (Q-k^)+ a^a2(a2-ai)9 V^k^(Q-2k^)]+ a^a| 9 

+ a^agk^U^^ 3 -9^7^[Q(2-5a^ag)-2k^(l-2aj^ag)] 

-7^ V^+ ■9^V^[Q(a^ag+2)+4k^(a^a2lag-a j^+7-12a^ag)-l)] 
+ 2-i)k^Bt2aj^cig9 ^Q+7^( l-a^ag)) + -v)^Q(a^a2Q+2aja2k2 

. (3a^-a2))l] + kV {Sa^ag { V^['?k2B(2Q-3k2)-2vV 

- k(a2-a^)J+ 'P^Q ^4V^ [k^(ag-a j^-2a^a2)+Q( (og- a^) 

. (oCi«2"^^"“i“2^^" ^“i®2'^^“2"“i^} +■^^^^{^(“2“®!^ 

. (3Q a a +2(k^-2Q))+29k^BQ(Q-k^(2a^ag+l)) 

J. ^ 

+ Y^[Q^+4k^a^(l-2ag)-k^Q(2(a^a2+ag-aj^)+5)]j^ - 
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- 3 f ( 2a g- 1 ) - S) B i 2-5a 2 ) + V ^ { a Q 

. (ag-a [Q(5a ^a^-2a\al-2}-^2k^ ( ioc^-a^) { l-5a ^ag ) 

+ a^a2(2a^a2-7)+l)]}]] + k^^-^-Sa^a^k V ( -p Q(,l-2a^a2) 

+ OgV^j+v'^ {■9^k[k^QUa^(i2(.3a2+a^-a^a2)-2) 

- a^CTgCag-a (,Q^+ 2 k‘^)]+V^-j^ k‘^7 +k' (2a ^ag (. 2 a 2 +a ^)- 1) 

- -9k^B[Q-2k^+2aj^ag(3k^-2Q)]+ 2 9 ^k^a ^[2k^Q(2a ^(Ba^-a 
-l)+4aj^a2k^+aga(l-Q)-4k^(l+4a|)]]]] - k^7® (og-a ^-2a ^Og) 

ikU ^-a^a2-?^®U®(ag-a^)-aj^agk®U^{ 

. [Q(2-a^ct2)-2k2]+ 9 V {k^[2(ag-ap (3a^a2-2) 

+ 2(2-90^02)]+ 2q(a^a2-l)] + 2 ol { 2 a ^- a ^) 

- Q(a2-a^)]- 2aj^a29^k -V^(a2-a^)-2k3B 9 [2a^ag 9 ^Q+V^]] 

- k^U^-[aj^ag{ 2kV(ag-a^)+7^ f 9^ [(Og-a^) (Sk^ (Q-k2)_Q2) 

+ 2k^ (4a ^ttg- 1 ) ( 2 ( q+k^ (k^- 1) ) +k^q)] -k- 9 k^B (og-a 

+ 9^7^k[2k(2a^ag-l)-3aj^a2q-9k^Bq(q-k^}(a2-a^)]} +V^ 

. 9^1^(3q-4k2)+72(a2-a^) [vV(2a^a2-l) + 9 k®B 

. ( 2 - 0 ^ 02 )+ 9^^[q(5a^a2-2)+4k^(l+2aj^a2(l~2a^ag)) 

+ 2(q+2a^ag)]]^j + k^U^ ^ 2a ^agV^^ [9 ^ (a ^ag-ag+a + 
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+ a^V )]- V k V 

^-;k2B^Q-2k2)U2-a^)]+.}2^k[a^a2iQ-k2)t,i_2ic2) 

- k^Qt l-2a ^a^)i2{a^-a ^)+a ^a^)]+2v\^ 

. [ ( 2k 2+ Q ) ( k 2-a ^a 2 ( 2k2+ Q )). 2a ^k^ J ^k^ [2 (i_4a ^ 

+ a^a2(5k2-4)) + Q(l-3a^a2)]]||| _ k^V® ( l-2a ^a^)] 

- a ^agS) k®U® ^ 2a -9 ^Q+V^ ( 2-5a 1*^2 L^ ^ ^ ) 

+ 2k^B-?+372]]' + k%®{_^kV(2-7a^a2-^2a|a|)-a^a2-9 

‘ i “ 1^2^ ^“ 2 ““ 1^ ^ ) +2 9 k^B i Q ( 1-a ^Og ) -k ^ ) 

+ '9^ Q^(2-5a^a2)+2Qk^ [BCug-a^) Ca ^^ag-D-ea^ag-l] 

- 8k^a^]] - Sa^Og^^ltQCag-a^)} + a^ag^k^BV^ 

+ ■v'^Q[V^(3a^a2Q-4k^(a2-a^) ( l-2a^ag) )-2a ^agk(ag-a ^) 

+ Ba^agT? Qlc^(Bk-2-? )] +9 V^[kta2-a^)(a^ag-2)+v2(2k2a^ 

. (2-3a^)+Q)+29k2BQ(l-a^a2)]]j +kV | V^k^BCag-cr ^ 

+ aj^ag(Q^+2k'^)]+kV^(a2-ai) l2ii^{l-2a^a^)+a^a^{a^-a^) 

. (Q-Dl+a^ttgkQ -i)^(a2-a^)[2Q(l-a^a2)+k^(3a^a2-2)] 

+ { k^y [(Og-a ^) (a ^ag- l)-5a j^ttg] +2k'^ [(Og-a ^) ( l-5a 3 ^ot 2 ■*■ 
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+ a^aglEk^+l)] + 2k^ S> | a [4a l+Q(i-3 (a^-a j) 

- 6a^a2)]+ 

+ Qi l-2a^a2(.3-4a^a2) ) + (,a2-aj_) Ua^ag 
-r k‘^BY^+2->)k7^[2a^a2(2k^+g) + ta2-a^)(.k^+Q)] 

- a^a29 [(Q-i^^)U2-a^)-tQ-2k2}]j] 
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Abstract 

In this paper, the VLSI complexity for the problem of merging two sorted arrays is studied. Let 
A and T be the chip area and computation time of a VLSI circuit that merges two sorted arrays 
of length V, by deriving the information transfer between the two sides of a partition of the 
chip, It is shown that log^ N) This result indicates that the inherent computational 

complexities for sorting and merging are asymptotically identical m the VLSI mode, unlike in 
the conventional model of sequential computation 

(Keywords * VLSI algorithms/parallel mergmg/area-time tradeoffs/computational complexities) 

Introduction 

With the advent of Veiy Large Scale Integration (VLSI), the chip area and 
computation time have been reduced dramatically. As a result, new methodologies 
and techniques are employed to design and manufacture small and expensive chips. 
Since in VLSI the wires, as opposed to the processing elements, dominate the 
manufacturing cost and computation time, the chip area rather than the number of 
logical elements is counted as the main measure of the complexity ,of the computation. 
Details on VLSI technology are discussed in Mead and Conway’ . 

One important and interesting area in VLSI research is to obtain lower bounds for 
the area-time tradeoffs for certain problems. During the past few years, various 
models for VLSI computation were developed'■^ Based on these models lower 
bounds for area - time tradeoffs are obtained for problems such as discrete Fourier 
transform^, integer multiplication^’* and matrix multiplication and some related 
problems’, sorting* transitive function* and a number of predicates*’®’®’*”. 
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In this paper, the problem of merging two sorted arrays is considered. Following 
the VLSI models discussed in literature'’^’'*-". We derive the lower bound for the area - 
time tradeoffs for this problem. Let A and T denote the chip area and computation 
time required to merge two sorted arrays of length N, we will show that AT^ = Q 
log^ N)andAT^=Q. V), for 0 < a: < 1 . 

The Merging Problem and VLSI Model 

Let I>={Q,\,-->Kr-\} be a set of K elements. Each element in D is called a word 
and can be represented with rflog bits. Let. £)//= {(jJi.JCa, ^n) 1 eD} be 

the set containing all the sorted arrays of length N. Let |y4|| denote the number of 
elements contained in a set^, then it is well known combinatorial result that 



N . 


( 1 ) 


The merging problem can be characterized by the function F •. Dn ><■ Dm £> 2 v. 
such that F (xi,...^M.yu-, Tw) = {zu-^ 2 n\ where Z= (xu-^n), Y= (yu-d>N) Dm, Z = 
(Zu...,ZzM ) e D2N, and {xx,...,Xm} u {yi,...yw,} = {Zx,...,Z-iM}. Arrays AT and F are the 
two input arrays to be merged, and array Z is the sorted array resulting from the 
merge. 

The VLSI model used in this paper is adopted from those proposed in 
literature'-^’'*. A VLSI circuit in our model can be viewed as the layout of a graph in 
which nodes correspond to logical gates and edges correspond to wires. While the 
nodes serve as either processing elements or I/O ports, the wires are used for 
interconnection between gates. We assume that chip layout of a VLSI circuit is a 
convex planar region. A wire in layout has minimal width A, > 0 and can transmit 1 bit 
of information . It is assumed that at most v > 2 wires can overlap at any point in the 
layout of a VLSI circuit, i.e., there are at most v layers in a chip. Each node can store a 
bit and occupies at least p > area. As for the processing time, it is assumed that a 
minimal delay t is required to transmit one bit through a node; this includes the delay 
over the wires. We also assume that each input is read exactly once and read in form a 
distinct input port, and that the times and locations at which the inputs are supplied 
and output are delivered are fixed and independent of the input values. The above 
assumptions reflect the state of nMOS circuit design methodology and fabrication 
technology'. Details on these assumptions can be found in literature'-^’". 
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Lower Bounds 


Based on the VLSI model described in the previous section, we will derive in this 
section a lower bound for AT^, where A and T are the chip area and computation time 
of a VLSI circuit that merges two sorted arrays. 

Our approach is now a classical one and consists of partitioning the circuit into 
two regions with roughly equal sizes, and then derive the minima! amount of 
information necessary to flow between the two sides of the partition during the course 
of computation. This techniques was originally developed in literature'* and was used 
to obtain the area time complexity of various problems in literature^ '*'®. Consider a 
function F (xi_..., Xfj) = and a VLSI circuit implementing F. Let L and R be 

the two regions obtained by partitioning the circuit with a chord perpendicular to the 
diameter of the layout. We will denote this partition by (I, R). Since it is always 
possible to slide the chord appropriately such that the number of input variables in R is 
equal to 5 = LiV/2j, we will assume throughout this paper that this is true for any 
partition (L.R). Also, for the sake of clarity in the remainder of this section, we assume 
that each input variable has all its bits input on the same side of partition. To be 
general we would have to define that a variable is input on the side of the partition 
where the majority of its bits are input. Since we would still consider at least half of 
the total number of bits, this would only change our result by a factor of 2. 

LetXfi = (jc,i,...,x„) ij <. . .< is, be the variables of Zin i? and Zr = (zji,...,Zjr), be the 
variables of Z in R, then the information flow from X to i? is defined in the following 
definition. 


Definition 1 : Let 




( 2 , Z,)| 


— C,,...,X„ — Cj, 


where Cr(ci,...,c^) e A be the restriction of F to Zr obtained by setting the input 
variables to the constants Ci,...,Cs respectively. The information -flow from X 

to R is defined by : 

I (F) = 

Definition 2 : The cross-flow I {F) that measures the amount of information that must 
flow between X and R during the computation of F is defined by 
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^L-*R ~ 


log max 


F =-C 

^ Xr 
CeDs 


( 2 ) 


The information-flow from R to L, ^r-^l is defined similarly. 

The cross flow 1(F) is related to ^4 7^ by the following lemma originally due to 
Thompson'*. 

Lemma 1 : The are A and the computation time T of a VLSI layout satisfy : 


Ar> 








( 3 ) 


where "k, x and v are the parameters discussed in the previous section. 

Lets consider now the problem of merging two sorted arrays. Let F 

yfi) = (zu...^ 2 n) be the function that implements merging. For 1 < /« <^ N-l by 
choosing yi =^ 2 = . . .=>’m = 0 and ym+i = ... = J^T-l, the restriction of F, 


Fm := F-. (xi,..xw,0,..,0,ii:-l K-\) ^ (0,...,0,x,,...x;,, is:-l,..,is:-l), (4) 


is the function that shifts each . input Xj by m position. Note that there are 
m O's and N<-m(K-iys in both parenthesis. Moreover, the sub function 
id : (xu...^N)-^(zm+u-,Zm+M) is the identify function in that 2 „+i = Xj. The function Fq 
and Ff/ are defined similarly by setting all the y/s to (^T-l) and 0 respectively. The 
cross flow of id is derived in lemma 2. The following notation is useful in our 
derivations: Given a set of subscripted variables W {u,u...,u,p} we let S(U) denote the 
set of subscripts {ii,...,ip}. 

Lemma 2 ; Let (L,R) be a portion of the circuit that implements the identify 
function id (xi,...^)^) = (zi,...,Za^), where X= (xx,...^)^),, Z= (zi,...^}^) e D^s and z, = xy, 
for i^j<N. LetZi andZ/j (Xr u Z« =Z) denote the sets of variables of Z that are in L 
and R respectively. Similarly, let Zl and Zr and (Zi u Zr = Z) denote the sets of 
variables of Z that are in I and i? respectively. Let \\S(Xi) n S’CZi?)!! + \\S(Xr) n ,S(Zi)||, 
the cross-flow of id satisfies. 
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I (id) = Cl (k \og k), (5) 

If A" = for some fixed value 8 > 0. 

Proof : Assume S{Xi) n where < ... <ip. By grouping all the 

variables with consecutive indices together, the set {x,x,...^,p} can be partitioned into t 

disjoint sets : 

^ }= k •■■■’% k k-. k „ 



where pj is the number of the variables in the y'* group. Note that pi+ ...+p/ = p and 
ip\=i\^P\-\, ip^^p^ = k^,+p2-l,...,etc. 

Now lets consider By choosing proper constant values for the input 

variables whose indices are not in we can derive the lower bound for ^l-^r 

by counting the number of joints in the range of 



since by the definition of ^l-^r. 


h^R > log 2 }| Xj = CjJ g ip } 


The constant value Cj assigned to XjJ € is defined as follows : 

(i) Cj = 0, ify < iu 

(ii) Cj= a„, ifj ^ h, if /p,+ ^p„ <J<ip,^ 1 ^ ^ 


(iii) Cj = Oi, ifj > , 
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where 

-0 and K, = \{K-\)pji\ > 1. 

Note that the choice of the sequence {«„} is legal, since {a„} is increasing and 

laa lasr P 

with Xj,j g being fixed to constant values as above, the cardinality of the 

range of (z,i,...,z,p) \xj = Cj j g is equal to the number of solution to the 

following inequality : 


0<x<...<x, <a,<x, <..jc, <a, ^...<<3, , <x <...<x <a,, 

'1 >p[ 1 'pi+i >PI+P2 2 t \ +ft-t+l 'n+ +;v ' 


( 6 ) 


And the log of the number of solutions is a lower bound for 
From the fact that the inequality 


0 ^ Xi ^...< Xp < ^-1 has 


k + p-l'] 


\ F J 


solution, the total number of solutions to (9) is 



%+ft-r 


X,+A-f 

^ A / 

k A ^ 


^ A y 


thus we have 




(K 


j-^Pj 

Pj 


■I') 


V 
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= ^ log 


^ S log 


fe] 

UJ 


pj 


^LPj log 

isja p 


= plog 


K-l 

p 


= a ip log p). 


The last equality follows from the fact that p<N and the assumption that K = 
for some fixed value s > 0. 

Similarly if ||>S'(Arjj) n S(Zi)\\ = q (q=^k-p), then we can show 


^R^L=Q.(q log q). 


(7) 


Since X^)=jc log x is a convex function, we have 


plogp + qlogq ^ p + q , p + g .g. 

2 2^2 

This gives I (id) = Q(A: log k). (9) 

Consider now the merge problem. Let (L,R) be a partition of the circuit that 
implements F(X,Y) = (Z) where X, Y e Dn aad Z e Dw We will introduced some 
useful notations as follows. Let F„ be the function that shifts each input variable Xj by 
m positions as defined in (6), and X^ X^ Zi and Z« be defined the same as in lemma 2. 
If {x,} / 6 / is a set of subscripted variables, then defined = {X,^m} i e I- Let 'F (F„) 
be defined by 
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Note that T (F„) measures the number of variables that "cross" through the 
partition of the circuit, i.e. the number of variables whose input ports and output ports 
lie in different sides of the partition. 

The following lemma shows that there exists at least an m., \ ^ m < N, for which 
the number of variables crossing the partition is Q(A0. 

Lemma 3 : There exists at least &nm, \ <m<N, such that 'P (F^) > LaVSJ. 

Proof : Let's consider the [^21 output variables xi_w/ 2 >i,. .,.xn and the [^21 output 
variables 

Only for each such Zj and x„ Xj is mapped to Zj by Fn where 1 < m = j-\ ^ N. 
Therefore, by counting through all the Fm, \^m<N, we have 


X'P(Fjs(rF/2l)min(/,/-), (11) 

X'F(Fjs(riV/2l) \NI2\I2\ (12) 

ISm^N 


and there must exist snm, \ <m<N, such that T (F») > LaVSJ . 

where / and r are the number of input variables x, LM2J +\<i<N, which are in L and 
R respectively. Since LM2J ^l,r <\nI2\, we have 

E'P(Fm) ^ 1 < 7w ^LM2J LF/ 2 J/ 2 J, and there must exist an m, \<.M<N such that 

'P(F„)>LM8J 

Combining lemma 2 and lemma 3, it follows that 

/(F„)0(FlogF). (13) 

Theorem : If /I is the chip area and T is the computation time required to merge 
two sorted arrays of lengths N, then 

(14) 


^7^ == Q {N^ log N) and 
= Cl log^* N), 


( 15 ) 
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for 0 < jc < 1 , if ^ for a fixed value e > 0. 

Proof : Since Fean realize each Fm,l <m<N, equation (14) follows directly from 
lemma 1 and equation (15). 

As for equation (15), consider the time when all the rA72l input variables 
have been read. None of the [N/2 \ variables Zw+i,...,ZwTa 721 have been 
output, since each x, can be routed to each Zj by some F^. These rM2l output variables 
can assume 2^^^^ distinct values; the circuit is therefore required to memorize I' N/2^, 
bits of information in order to assume 2^^^^ internal states. From the assumption in 
section 2 that each node occupies at least p area, we conclude that 

A>prM2l. (16) 

Equation (15) then follows directly from equation (14) and (16). 

Conclusion 

Following the VLSI model described in section 2, we have derived in the previous 
section lower bounds for the area-time tradeoffs for the merging two sorted arrays of 
lengths N is AT^ = Q log ^ A), For 0 < x < 1, F = for a fixed value s > 0. 
This shows that the inherent computational complexities for sorting and merging are 
asymptotically identical in the VLSI model'^ unlike in the conventional model of 
sequential computation. 

Although we required, in the proof given in the paper, that such input appear on a 
distinct input port, this assumption can easily be relaxed and it can be assumed that 
I/O ports are multiplexed. The proof would then follow arguments similar to the ones 
given in’^ for examples. 

Consider the general case of the merging problem in which the lengths of the two 
arrays to be sorted may not be the same. Lets assume that the lengths of these two 
arrays be M and N, where M'^N. From our results, it can be shown that the lower 
bound on area-time complexity for this general problem is ^7^ = f2(iV*^* log^ N) for 
0 < X ^ 1 . This lower bound is however not tight for the case when M is very large, 
compared with N. A t 5 q)ical example of this is the case when A=l. In that case, the 
merging problem is equivalent to the problem of inserting an element into a sorted 
arrays. An interesting open problem is therefore to obtain better lower bound for this 
general case of merging problem. 
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Equilibrium studies on mixed ligand complexes- 
Evidence of stacking interactions-solution studies 

Received March 14, 2002, Revised September 6 , 2002, Accepted February 4 , 2003 

Abstract 

Interaction of bivalent metals with maltol or kojic acid and selected ammo acids (Alanine, P- 
phenylalanine, L-tyrosine, L-DOPA and tryptophan) has been investigated by potentiomctnc 
pH measurements at 30 °C and p = 0 1 mol dm'^ (KNO 3 ). The stability order of these ternary 
complexes is in the order Ph alanine < tyrosine < DOPA < Tryptophan The extra stabilisation 
and above sequence of order is attributed to stacking interactions in ternary complexes 

(Keywords: potentiometry / ternary complexes / formation constants / maltol / kojic acid / 
amino acids). 

In continuation to our earlier work'‘^ on ternary metal complexes of y-pyrones, we 
now report the formation constants of ternary complexes of type MAL [where M = 
Co(II), Ni(n), Cu(n) and Zn(n); L = 3-hydroxy-2-methyl-pyran-4-one (maltol) or 5- 
hydroxy-2(hydroxymethyl)pyran-4-one (kojic acid); A = amino acids {viz, alanine, 
phenyl alanine, tyrosine, tryptophan and 3, 4 dihydroxyphenyl alanine (DOPA)] 
which have been determined potentiometrically at 30°C and p = 0.1 mol dm'^ (KNO^) 
in aqueous medium. 

The ligands, maltol and kojic acid are obtained from Lancaster (U.K). The amino 
acids are obtained from Sigma Chemical Company. The metal salts of Co(II), Ni(II), 
Cu(n), and Zn(II) are of AnalaR grade and their solutions are standardized 
volumetrically by titration with the sodium salt of EDTA.'' Potentiometric 
measurements are carried out at 30®C in aqueous medium (with 0.1 M KNO^ as 
background electrolyte) using control dynamics pH meter. The pH meter readings are 
plotted against a (moles of base added per mole of ligand) or m (moles of base added 
per mole of metal ion). Calculations were carried out with the help of computer 
program - BEST.^ The calculated ternary stability constants are reported in Table 1. 
Table 2 lists the AlogAT values along with AAlogA^ values. The parameters, Alog/C and 
AAlog^ are used to characterize the stability of ternary complexes formed. The extra 
stability found in ternary complexes is due to stacking interactions. 

The stacking interactions, can be quantised with suitable parameters associated 
with the following equilibrium: 
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Open form Closed form 

Table 1- Logarithms of stability constants (log ternary complexes 

[T= 30®C, p = 0.1 mol dm'^ (KNOj), aqueous medium] 


Secondary 

Ligands (A) 

Co(n) 

Ni(II) 

Cu(II) 

Zn(n) 


Kojic acid as primary ligand (L) 


alanine 

8.65 

10 06 

14.05 

9.56 

Ph.alanine 

8 53 

9.78 

13 78 

921 

Tyrosine 

8 02 

9.43 

13 15 

9 08 

DOPA 

8.65 

9 97 

14 16 

9.23 

Tryptophan 

8.70 

9.99 

13 90 

9 58 


Maltol as primary ligand (L) 


alanine 

9.28 

10 77 

15.14 

10.26 

Ph alanine 

9.20 

10.55 

14 91 

9.99 

Tyrosine 

8 70 

10 20 

14.26 

9 83 

DOPA 

9.34 

10.75 

15,31 

10.01 

Tryptophan 

9.35 

10.71 

15.04 

10.32 


Sigma fit values are m the range i 0 05 to 0.08 log units 
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The concentration independent and dimensionless constant K for intramolecular 
equilibrium is calculated by the relationship: ’ 

1 (i) 

The parameter AAlogA!^ is given by relationship 

AAlogA: = (Alog/Tj - AlogA:^) (2) 

2 + 

where AIogiiT^ is for the reaction [M + L + Ala] and AlogA:^ relates to the reaction 

2 + 

[M + L + other amino acids]. The values thus obtained are given in Table 2. 

Table 2-Companson of AlogK and AAlog^ values for ternary complexes 
[r= 30®C, p = 0 I mol dm"^ (KNO 3 ), aqueous medium] 


Secondary 

Ligands (A) 

Co(n) 

Ni(II) 

Cu(II) 

Zn(n) 


Kojic acid as primary ligand (L) 


alanine 

-0.38 

-0 37 

-0 45 

-0 39 

Ph alanine 

-0.32 

-0 27 

-0 31 

-0.33 


(+0.06) 

(+0 10 ) 

(+0 14) 

(+0 06) 

Tyrosine 

-0.29 

-0.25 

-0 26 

-0 26 


(+0 09) 

(+0 12 ) 

(+0.19) 

(+0.13) 

DOPA 

-0.26 

-0.16 

-0 21 

-0 26 


(+ 0 . 12 ) 

(+ 0 . 21 ) 

(+0.24) 

(+0.13) 


Table 2 Contd 
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Table 2 Contd 


Tryptophan 

-0 20 

-0 12 

-0 15 

-0 18 


(+0 18) 

(+0 25) 

(+0 30) 

(+0 21) 



Maltol as primary ligand (L) 



alanine 

-0 35 

-0.30 

-0 40 

-0 32 

Ph alanine 

-0 25 

-0 14 

-0 22 

-0 18 


(+0 10) 

(+0 16) 

(+0 18) 

(+0.14) 

Tyrosine 

-0.21 

-0 12 

-0.19 

-0 14 


(+0 14) 

(+0.18) 

^ (+0 21) 

(+0 18) 

DOPA 

-0.17 

-0 02 

-0.10 

-0.11 


(+0 18) 

(+0 28) 

(+0 30) 

(+0 21) 

Tryptophan 

-0.15 

-0 04 

-0 05 

-0.07 


(+0.20) 

(+0.26) 

(+0.35) 

(+0 25) 


Alog K = [log + log^j^ J. The AAlogA^ values are presented in paranthesis. 

AAlogA^ = [AlogATj - Alog/C 2 ], where Alog is for the reaction ( M + L + ammo acids except alanine) 
and Alog /<^2 for the reaction (M + L + alanine). 

Alanine is taken as a reference for zero based stacking interactions 

The percentage of stacked isomers in the systems is also calculated with the 
relationship; 

%(MLA)s, = [100 * / (1 + (3) 

The values of Ki and % (MLA)st for these systems are presented in Table 3. The 
perusal of parameters associated with stacking interactions reveal the following 
observations : 
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Table 3-Parameter showing the extent of intramolecular stacking interaction in temar>' complexes 
[T= 30®C, |j. = 0 1 mol dm'^ (KNO3), aqueous medium] 


Secondary 
Ligands (A) 

Parameters 

Co{II) 

Ni(II) 

Cu(ll) 

Zn(II 



Kojic acid as primary ligand (L) 



Ph alanine 

K, 

0 148 

0 258 

0 380 

0 148 


%(MLA)st 

12.89 

20 50 

27 54 

12 89 

Tyrosine 

K. 

0 230 

03182 

0 548 

0 348 


%(MLA)st 

18 70 

24 13 

35 44 

25 82 

DOPA 

Kj 

0318 

0 622 

0 738 

0 479 


%(MLA)st 

24 13 

38 35 

42 46 

32 38 

Tryptophan 

K, 

0513 

0 778 

0 995 

0 622 


%(MLA)st 

33.91 

43 75 

49 88 

38 35 



Maltol as primary ligand (L) 



Ph alanine 


0 258 

0 445 

0.513 

0 380 


%(MLA)st 

20.50 

30.79 

33 91 

27 54 

Tyrosine 

K. 

0.380 

0.513 

0.622 

0513 


%(MLA)st 

27 54 

33 91 

38.35 

33 91 

DOPA 


0 513 

0 9054 

0 9952 

0 622 


%(MLA)st 

33.91 

47.51 

49.88 

38 5 

Tryptophan 


0 585 

0 8197 

1 238 

0.778 


%(MLA)st 

36 91 

45 04 

55 32 

43.75 


^Benzimidazole as primary ligand (L) 


Ph alanine 

Kj 

1.137 

1.691 

1.344 

1 041 


%(MLA)st 

53.20 

62 83 

57.33 

51 00 

Tyrosine 


1.344 

1818 

1 884 

1 238 


%(MLA)st 

57.33 

64 51 

65 33 

55 31 

DOPA 

K, 

1 187 

1691 

1 454 

1089 


%(MLA)st 

54 27 

62 83 

59.25 

52.13 

Tryptophan 


1.454 

2 235 

2 235 

1 630 


%(MLA)st 

59 25 

69 08 

69.08 

61.97 

* Data for benzimidazole is abstracted from the reference. Usha Rani, G & DD Manwal (1998) 


J.Electrochem Soc India, 47»4 • 224 
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(a) The stacking interactions decrease in the order: tryptophan > tyrosine > 
DOPA > ph.alanine. This order can be correlated with the size of the 
aromatic side chain Stacking in ternary complexes containing tryptophan 
results in maximum^ stabilization since the surface area of the indole ring is 
large and permits for efficient stacking with kojic acid or maltol. The phenyl 
rings of ph.alanine and tyrosine are small and hence lower efficiency of 
stacking. 

The ligand DOPA is an ambidentate ligand. DOPA can bind to metal ion either 
like glycine using amino carboxylic group or it can bind like pyrocatechol using 
both oxygens present on phenyl ring. It is well known that in lower pH, DOPA 
binds to metal ion through amino carboxylate epd and at higher pH it binds like 
catechol. However as observed by Bhattacharya, the present study also shows that 
DOPA binds in ternary complexes from the amino carboxylate ends over the pH range 
of 4-7. 


(b) The % of stacking interaction decreases in the order Cu(II) > Ni(II) > Zn(II) 

> Co(II) with respect to the metal ions. This order is in confirmity with 
Irving Williams natural order of stabilities*. 

(c) The stability sequence with respect to primary ligand is in the order: Maloti 

> Kojic acid. This sequence is in line with basicity order. Maltol is more 
basic than kojic acid hence forms more stable complex. 

(d) The present communication deals with stacking interaction between 
aromatic moiety of one ligand and n electron cloud of another ligand. 
Earlier we have published a paper in which two ligands have aromatic 
moieties (eg. benzimidazoles as primary ligand and amino acids as 
secondary ligand). We thought it worthwhile to compare the extent of 
stacking interactions in these two different systems. Hence, we have 
reproduced the data of benzimidazole in Table 3. Comparison of % of 
stacked form from table 3 reveals the order : [benz-Cu-tryptophan](69.08%) 

> [maltol-Cu-tryptopahn](55.32%) > [kojic acid-Cu-tryptophan] (49.88%). 
Complex containing benzimidazole forms much better stacked complex 
than ternary complex of maltol or kojic acid. The reason is simple, the 
benzimidazole is aromatic whereas maltol or kojic acid are non-aromatic. 

Mohd. Zakee thanks the U.G.C. and the Principal, Mumtaz College for awarding 
teacher fellowship under F.I.P. 
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